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Abstract 


^An  investigation  was  conducted  to  examine  the  merits 
of  an  estimation  technique  involving  linear  interpolation 
to  estimate  Kolmogorov-Smirnov  (K-S)  critical  values 
when  the  scale  and  location  parameters  of  the  hypothesized 
distribution  are  unknown.  The  purpose  of  the  linear  esti¬ 
mation  technique  is  to  reduce  the  number  of  Monte  Carlo 
generated  samples  necessary  to  produce  useful  critical 
values  for  the  K-S  goodness-of-fit  test. 

Also,  different  plotting  positions  were  studied  to 
ascertain  which  plotting  positions  used  in  calculating 
and  plotting  the  K-S  test  statistic  values  provided  the 
best  critical  value  estimation.  \ 


The  distribution  used  asjtne  hypothesized  distri¬ 
bution  was  the  Type  I,  extreme  value  distribution  (lar¬ 
gest  extreme  value). 

C-v 

In  addition,  a  power  study  was  performed  which  com¬ 
pared  the  power  of  the  true  critical  values  against  the 
power  of  the  estimated  critical  values. 

( i 

The  following  are  the  major  results . Useful  critical 
values  were  found  with  the  linear  estimation  technique 
using  relatively  few  Monte  Carlo  generated  samples. 


viii 


v;. 


•  4  r. 


Further,  the  plotting  positions  found  to  be  the  best 
in  calculating  the  K-S  test  statistic  values  and 
plotting  these  values  were  respectively  i/n  and 
(i-.5)/n  where  i  is  the  ith  ranked  point  in  a  sample 
of  size  n. 
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A  MONTE  CARLO  TECHNIQUE  USING 
LINEAR  INTERPOLATION  TO  GENERATE 
MODIFIED  KOLMOGOROV-SMIRNOV 
CRITICAL  VALUES  FOR  THE 
EXTREME  VALUE  DISTRIBUTION 


I.  Introduction 


The  extreme  value  distribution  has  applications  in 
many  fields;  however,  its  application  to  reliability 
(Ref  12),  particularly  failure  due  to  corrosion,  and 
meteorology  (Ref  11),  make  it  important  to  the  Air  Force. 
Because  aircraft  reliability  and  meteorological  phenomena 
is  of  vital  interest  to  the  Air  Force,  data  pertaining  to 
both  areas  is  collected  by  the  Air  Force.  It  is  usually 
useful  to  determine  the  underlying  distribution  of  the 
reliability  or  meteorological  data.  A  "goodness-of-fit" 
test  can  be  performed  to  determine  the  underlying  distri¬ 
bution  for  the  data. 

A  goodness-of-fit  test  is  essentially  a  comparison 
between  a  hypothesized  distribution  and  the  frequency 
distribution  of  the  sample  data.  If  the  sample  data  "fit 
the  hypothesized  distribution  adequately,  then  the  hypo¬ 
thesized  distribution  is  used  as  the  representative  for 
the  underlying  distribution. 

The  Kolmogorov  Goodness-of-fit  Test 

The  Kolmogorov  goodness-of-fit  test,  which  will  be 


described,  is  often  referred  to  as  the  Kolmogorov-Smirncv 


r 


one-sample  test.  In  this  work,  the  test  shall  be  referred 
to  merely  as  the  Kolmogorov-Smirnov  (K-S)  test. 

The  purpose  of  the  K-S  test  is  to  determine  if  a 
hypothesized,  completely  specified,  continuous  distribu¬ 
tion  F^(x),  is  the  underlying  distribution  of  a  set  of 
sample  data,  ( x^  ,  x^  ,  .  .  .  ,  xn )  .  This  is  carried  out  by  com¬ 
paring  the  cumulative  distribution  function  F^Cx)  of  the 
hypothesized  distribution  to  the  cumulative  step  function 
Sn(x)  =  i/n;  i  equals  the  number  of  data  points  less  than 
or  equal  to  x,  and  n  is  the  number  of  data  points  in  the 
sample.  If  Fh(x),  the  proportion  of  the  population  defined 
by  the  hypothesized  distribution  having  values  less  than  x, 
is  adequately  close  to  the  step  function  Sn(x),  then  the 
hypothesized  distribution  is  accepted  as  the  underlying 
distribution.  A  graphical  view  of  this  comparison  is 
presented  in  Figure  1. 


1.0 

.5 

Fh(x)  r1 

✓X  DMAXy 

X  S  (x) 

/  n 

X 

Fig  1.  Graphical  View  of  the  K-S  Test  (Ref  4:346) 
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One  determines  if  Sn(x)  and  F^(x)  are  adequately 

close  by  using  the  maximum  vertical  distance  between 

F,  (x)  and  S  (x)  as  a  test  statistic.  DMAX  represents 
n  n 

this  maximum  vertical  distance  in  Figure  1.  In  mathe¬ 
matical  terms,  the  maximum  vertical  distance  is 


DMAX 


MAX 


F,  (x)  -  S  (x) 
h  n 


If  the  data  used  to  calculate  DMAX  is  from  the 
hypothesized  distribution,  the  distribution  of  DMAX  is 
known  and  is  independent  of  Fh(x);  however,  the  distri¬ 
bution  is  dependent  on  n  (Ref  2).  For  various  values  of 
n,  certain  DMAX  values  are  found  in  standard  tables  of 
critical  values  (Ref  20:70).  If  the  DMAX  statistic  for 
a  sample  of  size  n  is  larger  than  the  appropriate 
critical  value,  the  hypothesized  distribution  is  rejected 
as  being  the  underlying  distribution. 


Problem 

The  standard  tables  of  critical  values  for  the  K-S 
test  can  only  be  used  when  Fh(x)  is  completely  specified. 
In  other  words,  all  parameters  associated  with  F^(x) 
have  known  values.  However,  the  probability  integral 
transformation  used  by  David  and  Johnson  (Ref  8)  allows 
critical  value  tables  for  the  X-S  test  to  be  made  for 
any  particular  distribution  when  the  scale  and  location 
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parameters  of  F^(x)  are  unknown  and  estimated  from  the 
sample  data.  The  distribution  of  DMAX,  or  any  test 
statistic  based  on  the  cumulative  distribution  function, 
is  independent  of  the  scale  and  location  parameters, 
but  it  is  dependent  on  the  functional  form  of  the  hypo¬ 
thesized  distribution,  Fh(x)  (Ref  8:190).  Thus,  when 
either  the  scale  or  location  parameters,  or  both,  must 
be  estimated,  the  standard  critical  value  tables  cannot 
be  used,  and  new  critical  values  must  be  found  for  each 
hypothesized  distribution.  Lilliefors  has  tabulated 
critical  values  for  the  normal  distribution  with  unknown 
means  and  variance  (Ref  17).  Lilliefors  has  also 
generated  tables  for  the  exponential  distribution  with 
the  mean  unknown  (Ref  l8).  Tables  of  critical  values 
have  also  been  produced  for  the  three-parameter  Weibull 
distribution  and  the  gamma  distribution  where  the  scale 
and  location  parameters  in  each  were  estimated  (Ref  7). 

The  cumulative  distribution  function  of  the  extreme 
value  distribution,  which  will  be  defined  fully  later, 
has  associated  with  it  only  two  parameters.  These  para¬ 
meters  are  scale  and  location  parameters.  Hence, 
critical  value  tables  can  be  generated  for  the  extreme 
value  distribution  with  unknown  parameter  values. 
However,  this  table  can  only  be  created  through  Monte 
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Carlo  simulation  (computer  generation  of  an  empirical 
distribution  of  DMAX  values)  which  requires  a  large 
amount  of  computer  time. 

The  large  amount  of  computer  time  is  necessary  to 
generate  a  sample  size  of  DMAX  values  large  enough  to 
create  an  empirical  distribution  which  will  yield  valid 
critical  values.  One  of  the  research  questions  of 
this  study  is:  Is  there  a  technique  whereby  the 
necessary  sample  size  of  DMAX  values  can  be  lessened? 

Further,  the  cumulative  step-function  Sn(x)  defined 
earlier  is  also  referred  to  as  a  plotting  position  or  a 
plotting  convention.  It  is  the  plotting  position  nor¬ 
mally  used  in  Eq  (1)  to  calculate  DMAX.  It  is  also  the 
plotting  position  usually  utilized  in  plotting  the 
DMAX  values  in  order  to  obtain  an  empirical  distribution 
for  the  DMAX  statistic.  A  second  research  question  is: 
Is  Sn(x)  the  best  possible  plotting  position  to  use  in 
a  technique  which  would  lessen  the  necessary  number  of 
DMAX  values  created  through  Monte  Carlo  simulation? 

Background 

Some  preliminary  work  has  been  done  that  may  aid  in 
answering  the  questions  put  forth  in  the  previous  two 
paragraphs . 
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Johnston  (Ref  l4)  ,  working  with  component  relia¬ 
bilities,  as  opposed  to  defining  critical  values  of  a 
statistic,  found  good  results  from  empirical  distri¬ 
butions  formed  from  small,  Monte  Carlo-generated 
samples,  linear  interpolation,  and  linear  extrapolation. 
In  other  words,  linear  interpolation  between  the  data 
points  vaguely  outlining  the  empirical  distribution 
and  extrapolation  to  estimate  endpoints  allowed  a 
small  number  of  data  points  to  be  used  to  create  a 
usable  representation  of  the  true  distribution  of  the 
component  reliabilities.  A  small  number  cf  data 
points  means  a  small  amount  of  computer  time. 

Also,  other  plotting  positions  have  been  proposed 
(Ref  1:145;  15:300-301;  and  16:548-551)  as  better 
estimates  of  cumulative  distribution  functions  than 
S  (x).  One  of  the  other  plotting  positions  may  yield 
a  S-S  table  of  critical  values,  created  by  the  estima¬ 
ting  technique  previously  discussed,  that  is  more 
powerful  than  a  table  derived  using  Sr(x). 

Objectives 

This  work,  has  three  objectives.  The  first  objec¬ 
tive  is  to  discover  if  linear  interpolation  between  the 
Monte  Carlo-generated  data  points  (DMAX  values)  out¬ 
lining  the  DMAX  empirical  distribution  and  extrapolation 


used  to  obtain  endpoints  will  reduce  the  DMAX  sample 
size  required  to  find  suitable  critical  values. 

The  second  objective  is  to  investigate  the  possi¬ 
bility  that  plotting  positions  other  than  3n(x)  can 
produce  better  critical  value  estimates. 

The  third  and  final  objective,  which  is  a  product 
of  the  first  two  objectives,  is  to  fashion  a  table  of 
critical  values  for  the  estimating  technique  discussed. 

Overview 

This  study  provides  a  more  detailed  discussion  on 
several  of  the  points  brought  out  in  this  chapter.  Par¬ 
ticular  attention  is  given  to  the  description  and 
application  of  the  extreme  value  distribution  in 
Chapter  II. 

Chapter  III  gives  a  presentation  of  the  standard 
technique  for  finding  critical  values  along  with  a  des¬ 
cription  of  an  estimation  technique.  Also,  in  Chapter 
III,  plotting  positions  are  discussed. 

Chapter  IV  essentially  describes  those  activities 
undertaken  to  collect  and  analyze  the  data  used  in 
this  study. 

Chapter  V  tells  by  way  of  an  example  how  to  use 
critical  value  tables  for  the  K-S  test. 

Chapter  VI  presents  the  research  results. 


Chapter  VII  provides  concluding  remarks  about  this 
study  and  the  topic  area  it  covered.  Also,  comments 
concerning  how  this  study  could  have  been  improved  and 
how  this  study  could  be  expanded  are  offered  in  Chapter  VII. 

Finally,  the  appendices  present  the  reader  with 
tables  and  figures  which  are  useful  in  understanding  the 
development  and  results  of  this  thesis. 
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II.  The  Extreme  Value  Distribution 

Many  times  when  doing  a  statistical  analysis,  the 
major  interest  concerns  the  largest  or  smallest  extreme 
values.  Thus,  a  need  exists  to  understand  the  limiting 
form  of  the  distribution  of  the  smallest  or  largest 
extreme  value.  It  was  from  this  need  that  the  extreme 
value  distribution  arose. 

Evolution 

Consider  a  sample  consisting  of  n  data  points 
where  the  values  are  arranged  in  ascending  or  descending 
order.  Let  each  data  point  x  be  given  a  subscript  i  so 
that  i  indicates  order.  If  arranged  in  ascending  order, 
x^  is  the  smallest  value  and  xn  is  the  largest.  If 
several  ordered  samples  are  under  study  and  interest 
is  directed  to  one  of  the  extremes  in  each  sample,  then 
it  may  be  practical  to  have  some  knowledge  of  the 
probability  distribution  which  applies  to  such  a  set  of 
extremes . 

Consider  the  set  of  data  (x^ jXg, . • . »x  )  in  which  the 
elements  are  in  no  particular  order,  and  each  element  is 
defined  by  the  probability  density  function  g(x)  or  the 
cumulative  distribution  function  G(x).  Also,  let  W  be 
defined  as 
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(2) 


W  =  max{x. ,x.,....x  ) 

12  n 

Then,  the  cumulative  distribution  function  of  W,  in 
terms  of  x,  can  be  defined  as  H(x)  which  leads  to 

H(x)  =  P ( W <x ) 

=  (Xi£x,x2<  x . xn<x)  (3) 

Assuming  all  x^ ' s  are  independent,  Eq  (3)  becomes 

n 

H(x)  =  TTP(x.<x) 
i=l 

=  [G(x)]n  (4) 

The  probability  density  function  h(x)  is  then  the 
derivative  of  Eq  (4)  or 

h(x)  =  n[G(x)]n-1  g(x)  (5) 

Now,  consider 

Ln  H  ( x )  =  Ln  [  G  { x )  ] n 

=  nLn  G ( x )  ( 6 ) 

which  leads  to 

H ( x )  =  exp [ nLn  G ( x ) ] 

For  large  values  of  n,  the  exponent  of  e  in  Eq  (7) 
approaches  the  indefinite  form  infinity  times  zero. 
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Thus,  the  limiting  (or  asymptotic)  distribution  of  the 
largest  extreme  value  is  dependent  on  the  manner  in 
which  G(x)  approaches  unity  (Ref  23:66). 

Johnson  and  Kotz  (Ref  13:275-276),  explaining 
Gnedenko's  original  results,  show  the  correspondence 
between  the  behavior  of  G(x)  as  x  approaches  infinity 
and  the  Type  extreme  value  distribution  to  which  the 
limiting  distribution  belongs.  Three  Types  of  extreme 
value  distributions  exist. 

The  condition  which  leads  to  the  Type  I  extreme 
value  distribution  for  the  largest  extreme  value  is, 
as  x  approaches  infinity,  the  limit  of  ( 1  -  G(x)) 
converges  to  zero  as  quickly  as  the  limit  of  exp(-x) 
(Ref  23:66).  Several  parent  distributions,  the  distri¬ 
bution  defined  by  G(x),  meet  this  condition;  some  of 
them  are  the  normal,  log-normal,  exponential,  logistic, 
and  gamma  distributions  (Ref  9:5).  It  is  the  Type  I 
extreme  value  distribution  for  the  largest  extreme 
value  with  which  this  work  is  concerned. 

The  condition  leading  to  the  Type  I  extreme  value 
distribution  for  the  smallest  extreme  value  is,  as  x 
approaches  infinity,  the  limit  of  G(x)  converges  to 
zero  as  quickly  as  the  limit  of  exp(-x)  (Ref  15:^3) • 

To  give  a  broader  coverage  of  the  extreme  value 
distribution  it  will  be  noted  that  the  Type  II  limiting 
distribution  is  a  limiting  form  of  the  parent  distribu- 
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tion  of  the  Pareto  or  Cauchy  type.  The  Type  II  extreme 
value  distribution  has  not  been  applied  as  extensively 
as'  the  Type  I  or  Type  III  forms  and  is  considered  to 
be  of  lesser  importance  than  the  other  forms.  Further, 
the  Type  III  form  arises  when  the  range  of  the  density 
function  defining  the  parent  distribution  is  bounded 
from  below  (xlu),  and  the  cumulative  distribution 
function  of  the  parent  distribution  behaves  similarly 
to  (b-u)a  as  x  approaches  u  for  some  b>0  and  a>0 
(Ref  15:153)-  The  Type  III  extreme  value  distribution 
for  the  smallest  extreme  value  is  the  often  used 
Weibull  distribution,  which  has  been  applied  to  metal 
fatigue  and  fracture  problems. 

Description 

For  brevity  and  since  Type  I  is  the  one  most 
commonly  referred  to  as  the  "extreme  value  distribution," 
the  term  "extreme  value  distribution"  will  be,  from 
this  point  forward,  in  reference  to  the  Type  I  extreme 
value  distribution. 

The  cumulative  distribution  function  of  the  extreme 
value  distribution  (largest  extreme  value)  (Ref  13:272) 
is 


F(x)  =  exp[ -exp( -( x-u)/b ) ] 


(8) 


12 


and  in  the  same  manner  the  progression  was  made  from 
Eq  (4)  to  Eq  (5),  the  probability  density  function  is 
(Ref  22:9) 

f(x)  =  b_1exp[-(x-u)/b]exp[exp(-(x-u)/b) ]  (9) 

where  for  Eqs  (8)  and  (9)  -  Qo  <  X  <  ao  . 

The  cumulative  distribution  function  and  the 
probability  density  function  for  the  extreme  value 
distribution  (smallest  extreme  value)  are  respectively 
(Refs  1 5 «' 4-3 ;  22:9) 

F(x)  =  1  -  exp[-exp(x-u)/o]  (10) 

and 

f(x)  =  b-1exp[ (x-u)/b jexp[-exp(x-u)/b) ]  (11) 

where  for  Eqs  (10)  and  (11)  -o o<~^<oo. 

The  distributions  defined  by  Eqs  (8),  (9),  (10), 
and  (11)  have  a  scale  parameter  b  and  a  location  parameter 
u;  no  shape  parameter  is  involved.  The  graphic  represen¬ 
tation  of  Eqs  (9)  and  (11)  can  be  seen  in  Figure  2. 

Each  is  unimodal  with  its  mode  at  u  and  inflection 
points  at  (Ref  13:278) 

ui  bLn  [i(3+(5)*)l  (12) 
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Actually,  there  are  several  symmetry  relationships 
between  the  extreme  value  distribution  (largest  extreme 
value)  and  the  extreme  value  distribution  (smallest 
extreme  value).  These  are  seen  in  Table  I  and  Figure  2. 

As  can  be  seen  in  Figure  2,  both  distributions  are 
skewed  yet  mutually  symmetrical  about  u. 

Applications 

The  theoretical  basis  of  the  extreme  value  distri¬ 
bution  involves  consideration  of  the  parent  distribu¬ 
tion;  however,  application  of  the  extreme  value  distri¬ 
bution  does  not  always  demand  that  the  parent  distribu¬ 
tion  be  identified.  In  some  cases,  the  extreme  value 
distribution  has  been  applied  effectively  when  no  theor¬ 
etical  support  existed  for  its  use.  This  is  not  to  say 
that  it  should  be  applied  indiscriminately,  only  that 
it  is  not  always  necessary  to  analyze  the  nature  of  the 
parent  distribution.  This  has  allowed  the  extreme 
value  distribution  (largest  extreme  value)  to  be  utilized 
in  many  diverse  areas. 

Gumbel  (Ref  11)  has  applied  it  to  the  analysis  of 
structural  fatigue.  Gumbel  has  also  snown  applications 
of  the  extreme  value  distribution  (largest  extreme  value) 
to  flood  flows,  droughts,  and  other  extreme  meteorological 
phenomena  (Ref  12)  as  well  as  communication  systems 
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Table  I 


Characteristics  of  Extreme  Value  Distribution  (Ref  22:11) 


Characteristic 

Smallest  Value 

i 

Largest  Value 

Mean 

u  -  Sb 

Standard  Dev. 

b  Tf/-Je 

Median 

u  -  .3665 l  b 

Mode 

u 

S  =  .57721. . . 


Fig  2. 

Extreme  Value  Distributions 
of  Smallest  and  Largest  Values 
with  b  =  1,  u  =  0  (Ref  22:11) 


(Ref  10).  Bogdanoff  and  Schiff  (Ref  3 :  147-178)  have 
shewn  the  applicability  cf  the  distribution  to  the  pre¬ 
diction  of  the  force  of  earthquakes.  Johnson  and  Kotz 
(Ref  13:274)  have  referred  to  numerous  applications; 
microorganism  survival  times,  corrosion,  and  dielectric 
strengths  of  paper  capacitors.  Kapur  and  Lamberson 
cite  its  use  in  analyzing  gust  velocities,  gust  loads, 
and  landing  loads  for  aircraft  (Ref  15:154).  It  has 
been  used  by  Branger  (Ref  4:213-237)  to  estimate  the 
life  of  fighter  aircraft. 

Chapter  II  Summary 

A  need  to  analyze  the  extreme  values  of  a  sample 
was  the  impetus  in  the  evolution  of  the  extreme  value 
distribution.  In  the  early  theoretical  work  on  the 
extreme  value  distribution,  it  was  discovered  that  the 
extreme  value  distribution  has  three  forms.  These  forms 
are  called  Type  I,  Type  II,  and  Type  III.  The  Type  of 
the  extreme  value  distribution  is  dependent  on  the 
limiting  properties  of  the  parent  distribution. 

The  extreme  value  distribution  is  applied  widely. 

It  is  even  being  applied  successfully  in  circumstances 
where  the  theoretical  foundation  for  its  use  is  weak. 
Several  of  the  areas  in  which  the  extreme  value  distri¬ 
bution  is  applied  are  structural  fatigue,  flood  flows, 
and  gust  velocities. 
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III.  The  Techniques  and  Plotting  Positions 


The  standard  method  of  finding  critical  values  for 
the  K-S  test  when  the  hypothesized  distribution  is  fully 
specified  usually  involves  a  Monte  Carlo  technique 
which  generates  a  large  number  of  samples.  Normally, 

5,000  to  10,000  samples  are  generated.  From  each  of 
these  samples  a  DMAX  value  is  found,  thus  giving  a 
large  sample  of  DMAX  values.  With  such  a  large  number 
of  DMAX  values,  an  excellent  approximation  to  the  cumu¬ 
lative  distribution  function  of  the  DMAX  values  can  be 
made.  Actually,  with  so  many  values,  the  critical 
values  found  using  this  standard  technique  are  taken 
as  truth. 

On  the  other  hand,  the  linear  estimation  technique 
uses  a  drastically  reduced  number  of  samples  as  compared 
to  the  standard  technique. 

Additionally,  it  should  be  pointed  out  that  plotting 

positions  other  than  S  (x)  have  been  advanced.  The  use  of 

n 

an  alternative  plotting  position  may  offer  a  better 
linear  estimation  technique. 

Standard  Technique 

The  method  which  shall  be  discussed  is  the  method 
normally  used  to  obtain  critical  values  at  their  corres¬ 
ponding  significance  levels  for  the  K-S  test  when  the 
hypothesized  distribution  is  not. fully  specified. 
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For  a  fixed  sample  size  n,  a  sample  of  extreme  value 


(largest  extreme  value)  random  deviates  is  generated. 

The  scale  and  location  parameters  are  then  estimated 
since  the  hypothesized  distribution  is  not  fully  speci¬ 
fied.  DMAX,  which  was  introduced  in  Sq  (1),  is  calcu¬ 
lated  for  the  sample.  This  procedure  is  repeated  using 
a  different  set  of  sample  elements  each  time.  Assume 
the  procedure  is  repeated  5>000  times  thus  giving  a 
sample  consisting  of  55000  independent  DMAX  values. 

The  5,000  DMAX  values  are  ranked  and  the  80th,  85th, 

90th,  95tn,  and  99th  percentiles  are  found.  The  percen¬ 
tiles  are  interpreted  then  as  the  critical  values. 

Ranking  the  DMAX  values  and  finding  the  percentile 
values  is  analogous,  in  a  graphical  sense,  to  plotting 
the  DMAX  values  with  respect  to  the  vertical  axis  using 
Sn(x).  Thus,  the  vast  number  of  plotted  points  is  an 
approximation  to  the  plot  of  the  cumulative  distribution 
function  of  the  DMAX  values.  This  interpretation  can 
be  seen  in  Figure  3*  With  5»000  points,  the  curve  is, 
for  all  practical  purposes,  completely  defined. 

It  should  be  stated  that  associated  with  each 
critical  value  is  a  significance  level.  The  significance 
level,  often  denoted  merely  by  a,  is  the  risk  of  rejecting 
the  hypothesized  distribution  as  the  true  underlying 
distribution  when  in  fact,  it  is  the  underlying  distribu- 
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tion.  In  other  words,  the  DMAX  value  for  a  sample  of 
size  n  which  has  the  extreme  value  distribution  (largest 
extreme  value)  as  its  underlying  distribution  will  be 
greater  than  the  critical  value  (denoted  DMAXn,  a) 
associated  with  sample  size  n  and  significance  level  a 
with  probability  a.  Thus,  the  critical  values  corres¬ 
ponding  to  the  80th,  85th,  90th,  and  99th  percentiles 
have  significance  levels  which  are  respectively  .20,  .15, 
.10,  .05,  and  .01. 

The  Estimation  Technique 

The  technique  for  estimating  critical  values  is 
described  below.  It  is  detailed  in  terms  of  a  graphical 
procedure  for  reasons  of  clarity. 

Initially,  for  a  fixed  sample  of  size  n,  a  sample  of 
extreme  value  (largest  extreme  value)  random  deviates  is 
generated,  and  from  these  n  random  deviates,  the  scale 
and  location  parameters  are  estimated.  Using  these  esti¬ 
mated  parameter  values,  DMAX  is  calculated  for  the  sample. 
This  procedure  is  repeated  N  times  with  a  different  sample 
of  random  deviates  for  each  repetition.  The  N  DMAX  values 
are  plotted  with  respect  to  the  vertical  axis  using  a 
plotting  position,  and  lines  are  then  drawn  between  the 
plotted  points.  This  concept  can  be  seen  in  Figure  4. 
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s„(x) 


DMAX 


=  critical  value 


Fig  3.  Graphical  Representation  of  Standard  Technique 
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Fig  4.  Graphical  Representation  of 
Linear  Estimation  Technique 
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Further,  lines  are  extrapolated  outward  to  create 
endpoints;  this  can  be  seen  in  Figure  4  where  the  seg¬ 
mented  lines  represent  the  extrapolated  lines.  The  end¬ 
points  are  formed  at  T(DMAX)  =  1  and  T(DMAX)  =  0  where 
T(DMAX)  is  the  cumulative  distribution  function  of  DMAX. 

The  critical  value  DMAX  is  found  by  discerning 
what  DMAX  value  corresponds  with  the  vertical  axis  value 
1  -  a.  For  example,  the  critical  value  DMAXn  is 

found  in  Figure  4  by  deducing  the  DMAX  value  which  pairs 
with  F(DMAX)  =  .90.  All  critical  values  are  found  in 
this  manner.  As  can  be  seen  in  Figure  4,  the  bar-dot 
lines  which  "pick-off"  the  critical  value  intersect  not 
at  one  of  many  points,  as  in  Figure  3,  but  at  the  line 
drawn  between  two  of  merely  a  few  points;  thus,  the 
critical  value  is  estimated. 

Plotting  Positions 

In  the  description  of  the  standard  technique  for 
finding  critical  values,  it  was  alluded  that  the  plotting 
position  Sn(x)  is  used  as  an  approximation  to  F^(x)  and 
thus  is  used  in  the  calculation  of  DMAX.  And,  as  already 
mentioned  Sn(x)  is  used  to  "plot"  the  DMAX  values.  How¬ 
ever  $n(x)  is  not  the  only  plausible  plotting  position. 

An  alternative  may  prove  to  be  more  useful  in  the  linear 
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estimation  technique.  Blom,  Kimball,  and  others  have 
suggested  alternate  plotting  positions  with  which  to 
use  to  approximate  cumulative  distribution  functions 
(Refs  1:145;  15:301,  and  16:5^8-549). 

Five  alternative  plotting  positions  are  investigated 
in  this  work.  One  of  the  alternatives  investigated  is 


P(xi) 


i 

n+1 


(13) 


It  can  be  shown  that  Eq  (13)  is  the  mean  value  of  the  frac¬ 
tion  of  the  population,  defined  by  the  cumulative  distri¬ 
bution  function  F(x),  falling  prior  to  the  ith  ranked  data 
point  in  a  sample  of  size  n  ( Ref : 279-299 ) •  Eq  (13)  is 
often  called  the  mean  rank. 

Another  alternate  plotting  position  is 


F(x,  ) 

X 


i-«  3 

n+  .4 


(14) 


Eq  ( 1 4 ) ,  called  the  median  rank,  is  the  median  value  of 
the  fraction  of  the  population,  defined  by  F(x),  coming 


before  the  ith  ranked  data  point. 

The  other  plotting  positions  examined  are 


Chapter  III  Summary 

The  standard  technique  for  finding  critical  values 
requires  that  a  large  number  cf  DMAX  values  be  generated. 
Normally,  5,000  to  10,000  DMAX  values  are  used.  This 
large  number  of  DMAX  values  "plotted”  using  S_(x)  defines 

il 

the  DMAX  cumulative  distribution  so  closely  that  the 
critical  values  found  using  this  technique  are  considered 
true . 

The  linear  estimator  technique  uses  drastically 
fewer  DMAX  values.  These  DMAX  values  are  plotted  using  a 
plotting  position,  not  necessarily  Sn(x).  Lines  are 
drawn  between  the  plotted  DMAX  values,  thus  giving  a 
rough  outline  of  the  graph  cf  the  cumulative  distribution 
of  the  DMAX  values. 

Several  different  plotting  positions  have  been 
suggested  as  alternatives  to  S  (x) .  At  least  one  of  these, 
as  opposed  to  Sn(x),  may  supply  a  better  linear  estimating 
technique . 


IV.  Methodology 

The  following  sections  describe  the  overall  process 
that  was  undertaken  to  acquire  and  analyze  the  data  needed 
for  this  study.  A  detailed  description  is  given,  and 
then,  in  the  summary  section  a  step-by-step  overview  of 
the  methodology  process  used  in  this  work  is  presented. 

Initial  Process 

Initially,  the  standard  technique,  using  5 ,000  DMAX 
values,  was  employed  to  generate  critical  values  for 
n  =  10,  25,  and  40  at  the  significance  levels  .20,  .15, 
.10,  .05,  and  .01.  These  critical  values,  which  shall 
be  referred  to  as  the  standard  set  of  critical  values, 
were  taken  as  the  true  critical  values  and  used  for  com¬ 
parison  with  the  estimated  critical  values. 

Next  N  DMAX  values  were  generated  at  n  = 10,  25, 
and  40.  These  were  utilized  in  the  linear  estimation 
technique  to  estimate  the  true  critical  values  at  the 
significance  levels  .20,  .15,  .10,  .05,  and  .01.  The 
N  DMAX  values  were  calculated  using  Sn(x)  in  Eq  (1) 
and  plotted  using  Sn(x)  and  each  of  the  five  alternative 
plotting  positions  thus  giving  six  sets  of  estimated 
critical  values.  A  set  of  estimated  critical  values  was 
compared  to  the  standard  set  of  critical  values  by  cal¬ 
culating  the  mean  absolute  difference  between  the 
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elements  of  the  estimated  set  and  the  corresponding  ele¬ 
ments  of  the  standard  set.  An  additional  comparison  was 
made  for  each  estimated  set;  note  was  taken  of  the  lar¬ 
gest  absolute  difference  between  the  elements  of  the  es¬ 
timated  set  and  the  corresponding  elements  in  the  standard 
set.  This  comparison  procedure  involving  the  mean 
absolute  difference  and  the  largest  absolute  difference 
occured  for  N  =  100,  250,  500,  and  750. 

It  was  viewed  that  a  "good”  estimation  would  be 
approximately  a  mean  absolute  difference  of  .002  and  a 
largest  absolute  difference  of  .004.  It  should  be  stated 
that  the  two  values  given  as  measures  of  a  "good"  estima¬ 
tion  were  arbitrarily  selected.  It  was  also  felt  that 
going  to  N  larger  than  750  to  reach  the  measures  of  a 
"good"  estimation  would  signify  that  the  linear  estima¬ 
tion  technique  is  a  poor  estimating  technique.  This 
view  was  taken  since  the  primary  objective  of  the  tech¬ 
nique  is  to  markedly  reduce  N,  the  number  of  DMAX  values. 

From  comparing  each  of  the  six  sets  of  estimated 
critical  values  to  the  standard  set,  the  best  plotting 
position  with  which  to  plot  the  DMAX  values  was  chosen. 
This  plotting  position  was  then  used  to  plot  the  DMAX 
values  in  all  additional  implementations  of  the  linear 
estimation  technique. 
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Of  the  five  alternative  plotting  positions  to  3  (x), 
the  three  best  in  terms  of  plotting  the  DMAX  values  were 
chosen  to  be  analyzed  as  alternatives  to  S  (x)  in  the 
calculation  of  DMAX.  It  was  felt  that  the  plotting  posi¬ 
tions  which  best  estimated  the  cumulative  distribution 
function  of  DMAX  might  provide  the  best  estimation  of 
F,  (x),  in  this  case,  the  cumulative  distribution  function 

Ij. 

of  the  extreme  value  distribution  (largest  estreme  value). 
Thus,  three  new  sets  of  estimated  critical  values  were 
compared  to  the  standard  set  giving  a  mean  absolute 
difference  value  and  a  largest  absolute  difference  value 
for  each  estimated  set.  From  these  values  it  was  deter¬ 
mined  if  substituting  Sn(x)  with  some  other  plotting 
position  from  the  five  alternatives  set  forth  would  pro¬ 
vide  a  better  linear  estimation  technique. 

Essentially,  the  proceeding  process  gave  analysis  in 
three  important  areas  of  this  study.  First,  it  analyzed 
the  matter  of  whether  or  not  the  linear  estimation  tech¬ 
nique  is  a  "good"  technique  in  terms  of  providing  crit¬ 
ical  values  "close"  to  the  true  critical  values.  Secondly, 
it  investigated  other  plotting  positions  as  substitutes 
for  S  (x)  in  plotting  the  DMAX  values.  Finally,  it  studied 
the  possibility  that  at  least  one  of  the  alternative 
plotting  positions  should  replace  S  (x)  in  the  calculation 
of  DMAX. 
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One  further  topic  was  analyzed;  that  was  the  power 
of  the  linear  estimation  technique.  How  this  analysis 
was  done  is  discussed  in  the  following  section. 

Power  Study 

The  term  power  is  defined  as  the  probability  of 
rejecting  the  hypothesized  distribution  as  the  underlying 
distribution  of  a  sample  when  the  hypothesized  distri¬ 
bution  is  not  the  underlying  distribution.  From  this 
definition,  it  can  be  seen  that  a  goodness-of-f it  test 
with  high  power  values  is  desirable. 

Initially,  a  power  study  was  conducted  on  the 
linear  estimation  technique  (with  N  =  100,  250,  500, 
and  750)  formed  from  the  plotting  positions  that  were 
determined  to  give  the  best  technique. 

The  power  study  on  the  linear  estimation  technique 
was  executed  by  generating  2,500  samples  of  some  distri¬ 
bution  other  than  the  extreme  value  distribution  (largest 
extreme  value).  For  each  sample,  a  DMAX  value  was  cal¬ 
culated  using  the  plotting  position  which  was  chosen 
earlier  to  be  the  best  in  calculating  DMAX.  These  DMAX 
values  were  compared  to  the  appropriate  critical  values, 
and  the  number  of  rejections  of  the  extreme  value  distri¬ 
bution  (largest  extreme  value)  as  the  underlying  distri¬ 
bution  were  counted.  The  proportion  of  rejections  out 
of  the  2,500  trials  was  taken  as  the  power. 
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This  process  was  conducted  for  the  following  distri¬ 
butions:  the  two -parameter  Weibull  (shape  =  1,  scale  =  1), 

the  standard  normal,  the  two -parameter  log-normal  (mean  =  0, 
standard  deviation  =  1),  the  chi-square  (1  degree  of  free¬ 
dom),  and  the  chi-square  (4  degrees  of  freedom).  Power 
values  were  also  found  using  the  same  distributions, 
sample  sizes,  and  significance  levels  for  the  standard 
critical  values.  Using  Sn(x)  in  the  calculation  of  DMAX, 
2,500  DMAX  values  were  found  and  compared  with  the  appro¬ 
priate  standard  critical  values.  The  power  values  from 
the  estimated  technique  were  compared  to  the  power  values 
from  the  standard  technique. 

Generation  of  Random  Deviates 

In  this  study  it  was  necessary  to  generate  random 
deviates  from  six  different  distributions.  These  distri¬ 
butions  were  the  extreme  value  (largest  extreme  value), 
the  two-parameter  Weibull,  the  standard  normal,  the  two- 
parameter  log-normal,  the  chi-square  (1  degree  of  freedom), 
and  the  chi-square  (4  degrees  of  freedom). 

The  extreme  value  (largest  extreme  value)  random 
deviates  were  generated  with  the  use  of  two-parameter 
Weibull  random  deviates.  The  two-parameter  Weibull 
distribution  is  defined  by  the  cumulative  distribution 
function 

F(t)  =  1  -  exp[-(t/z)m]  ,  t  _>  0  (18) 
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where  the  shape  parameter  m  and  the  scale  parameter  z 
are  both  assumed  to  be  positive. 


If  the  transformation  x  =  -Ln(t)  is  made  when  t  has 
the  distribution  defined  by  Eq  (18),  then  x  is  distributed 
according  to  the  extreme  value  distribution  (largest 
extreme  value)  (Ref  5:2).  The  scale  and  location  parameters 
of  the  extreme  value  distribution  (largest  extreme  value) 
are  then  b  =  m-^  and  u  =  -Ln(z)  respectively. 

Thus,  in  order  to  generate  extreme  value  random 
deviates,  two-parameter  Weibull  random  deviates  were 
generated.  The  parameter  values  used  were  z  *  1  and  m  =  1. 
The  Weibull  random  deviates  were  generated  on  the  CDC  6600 
computer  by  the  International  Mathematical  and  Statistics 
Library  (IMSL)  subrouting  GGWIB.  The  transformation 
x  s  -Ln(t)  was  performed  on  the  two -parameter  Weibull 
random  deviates,  and  extreme  value  random  deviates  with 
u  =  0  and  b  =  1  were  then  produced. 

The  standard  normal  random  deviates  were  generated 
by  using  the  IMSL  subroutine  GGNML. 

The  two-parameter  log-normal  random  deviates  were 
generated  with  the  use  of  the  subroutine  GGNLG  which  is 
also  from  IMSL. 

The  chi-square  (1  degree  of  freedom)  random  deviates 
were  generated  through  the  use  of  standard  normal  random 
deviates.  It  can  be  shown  that  the  square  of  a  standard 
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normal  random  variable  has  a  chi-square  (1  degree  of 
freedom)  distribution  (Ref  21:217-218).  Using  that 
fact,  standard  normal  random  deviates  were  generated, 
and  these  random  deviates  were  squared.  Thus,  chi-squared 
(1  degree  of  freedom)  random  deviates  were  formed. 

The  chi-square  (4  degrees  of  freedom)  random  deviates 
were  also  generated  using  standard  normal  random  deviates. 
The  method  that  was  employed  is  the  same  as  that  used  by 
Littell,  McClave,  and  Offen  (Ref  19:265).  Each  chi-square 
(4  degrees  of  freedom)  random  deviate  was  constructed  by 
summing  the  squares  of  four  independent  standard  normal 
random  deviates. 


Estimation  of  Parameters 


The  scale  and  location  parameters  were  estimated  by 
using  recursive  equations  presented  by  Stephens  (Ref  24). 
These  equations  provide  maximum  likelihood  of  the  para¬ 
meters.  The  equation  for  the  location  parameter  is 


-1 


j  exp(-x./b) 
J 


u  =  -b  Ln  j^n 
and  the  equation  for  the  scale  parameter  is 


(19) 


b  =  n 


-1 


-j  XJ  " 


.x  .  exp  ( — x/b  )| 
’ J  J 


exp( -x/b ) 


-1 
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In  both  Eqs  ( 1 8 )  and  (IS),  n  is  the  number  of  elements 

in  the  sample,  and  x.  is  the  jth  element  in  the  sample. 

J 

Further,  all  summations  represent  sums  from  j  =  1  to  n. 

As  can  be  seen  in  the  equations,  b  must  be  solved  for 
recursively  to  some  arbitrary  tolerance  level  and  then  sub¬ 
stituted  into  Eq  ( 1 8 )  to  find  u.  The  tolerance  level  used 
in  finding  b  was  .00001.  That  is,  if  the  change  in  the 
previous  estimate  of  b  is  less  than  or  equal  to  .00001 
when  compared  to  the  current  estimate,  the  estimation 
process  is  discontinued,  and  the  current  estimate  is 
accepted.  The  maximum  number  of  iterations  allowed  was 
500;  however,  this  was  never  reached. 

It  should  be  noted  that  in  order  to  initiate  the 
recursive  process  to  find  b,  an  initial  guess  of  b  must  be 
made.  This  guess  was  b  =  1  in  all  circumstances. 

Computer  Programs 

The  computer  programs  utilized  in  this  work  can  be 
found  in  Appendix  G.  These  computer  programs  were  made  to 
be  adaptable  by  allowing  the  user  to  easily  change  program 
parameters.  In  addition,  the  programs  are  segmented  by 
comment-statement  headings  which  state  the  function  of  the 
segment  it  designates. 

Chapter  IV  Summary 

The  computer  programs  used  in  this  study  were  made 
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flexible  by  permitting  the  user  to  quickly  change  program 
parameters . 

Random  deviates  were  needed  from  several  different 
distributions.  The  distributions  were  the  extreme  value 
(largest  extreme  value),  the  standard  normal,  the  two - 
parameter  log-normal,  the  two-parameter  Weibull,  the  chi- 
square  (1  degree  of  freedom),  and  the  chi-square  (4  degrees 
of  freedom) . 

The  scale  and  location  parameters  were  estimated  by 
recursive  equations.  The  tolerance  level  used  was  .00001, 
and  the  maximum  number  of  iterations  allowed  was  500. 

The  other  actions  taken  in  the  process  of  acquiring 
and  analyzing  data  can  be  roughly  outlined  in  the  steps 
which  follow: 

1.  A  standard  set  of  critical  values  was  formed. 

2.  A  determination  was  made  as  to  whether  or  not  the 
linear  estimation  technique  was  a  "good"  estimation  tech¬ 
nique  . 

3.  A  determination  was  made  as  to  which  plotting 
position  best  plotted  the  DMAX  values. 

4.  Of  the  five  alternative  plotting  positions,  the 
three  best  in  terms  of  plotting  DMAX  values  were  analyzed 
as  substitutes  for  Sn(x)  in  calculating  DMAX. 

5.  It  was  determined  which,  if  any,  of  the  three 


plotting  positions  chosen  in  step  3  should  be  substituted 


for  S  (x)  in  DMAX  calculations. 

n 

6.  The  plotting  positions  selected  in  step  2  and 
step  4  were  designated  as  forming  from  the  plotting 
positions  studied  the  best  linear  estimation  technique. 

7.  A  power  study  was  performed  on  the  linear  estima¬ 
tion  technique  found  in  step  5- 


V.  Use  of  Critical  Value  Table 


The  process  required  to  use  a  K-S  critical  value 
table,  when  the  location  and  scale  parameters  are  estima¬ 
ted,  is  quite  simple.  The  DMAX  value  is  calculated  for 
the  n  data  points  which  make  up  the  sample.  This  DMAX 
value  is  compared  to  DMAXn  If  DMAX  for  the  sample  is 

greater  than  DMAXn  a,  the  hypothesized  distribution  is 
rejected  as  the  underlying  distribution  of  the  sample. 

A  better  explanation  can  be  provided  through  a  demonstra¬ 
tion  . 

An  Example 

The  example  which  follows  involves  the  use  of  the 
critical  value  table  generated  using  the  standard  tech¬ 
nique;  however,  the  critical  value  table  generated  with 
the  linear  estimation  technique  is  used  in  the  same  manner. 

Assume  the  sample  data  -.56,  1.54,  0.17,  3-42,  2.11, 
1.16,  4.Q2,  -1.23,  -1.71,  3-72,  4.23,  1.91,  3-37,  -2.36, 
and  2.13  has  been  collected.  Also  assume  that  the  extreme 
value  distribution  has  been  hypothesized  as  the  underlying 
distribution  of  the  sample  and  that  this  hypothesis  is  to 
be  tested. 

Before  the  DMAX  value  for  the  sample  can  be  found,  the 
scale  and  location  parameters  must  be  estimated  from  the 
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fifteen  sample  elements.  These  parameters  could  possibly 
be  estimated  by  hand  from  Eqs  (19)  and  (20)  but  only  with 
some  difficulty.  Fortunately,  Eikman  has  tabled  two  sets 
of  coefficients  which  can  be  used  to  calculate  unbiased 
nearly  best  linear  estimates  of  the  scale  and  location 
parameters  of  the  extreme  value  distribution  (largest 
extreme  value)  (Ref  9:32-423).  By  multiplying  each 
sample  element  by  the  appropriate  coefficient  and  summing 
the  products,  the  scale  or  location  parameter  can  be 
found  depending  on  which  set  of  coefficients  is  used. 

The  ease  of  computation  using  the  coefficients  more  than 
compensates  for  the  slight  sacrifice  in  accuracy  of  the 
estimate.  Now,  using  the  two  sets  of  coefficients  on 
the  sample  elements  reveals  estimates  of  2.14  for  the 
scale  parameter  and  0.34  for  the  location  parameter. 

Substituting  the  estimated  parameter  values  into 
the  cumulative  distribution  function  of  the  extreme  value 
distribution  (largest  extreme  value)  allows  Table  II  to 
be  constructed.  In  Table  II,  the  DMAX  value  is  promin¬ 
ently  indicated  by  asterisks.  From  Table  III,  the  stan¬ 
dard  table  of  critical  values,  it  can  be  seen  that  DMAX.,- 

15, .05 

=  .2194.  Since  the  DMAX  value  for  the  sample  is  less 

than  the  critical  value  DMAX.-  the  extreme  value 

15 , . U5 

distribution  (largest  extreme  value)  cannot  be  rejected 
as  the  underlying  distribution. 
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The  steps  necessary  to  use  a  K-S  critical  value 
table  for  the  extreme  value  distribution  (largest  extreme 
value)  when  the  location  and  scale  parameters  are  unknown 
are  these: 

1.  Estimate  the  location  and  scale  parameters  from 
the  sample  data. 

2.  Find  the  DMAX  value  for  the  sample. 

3.  Compare  the  DMAX  value  for  the  sample  to  the 
critical  value  DMAX 

n,a 

4.  If  the  sample  DMAX  value  is  greater  than  DMAX  , 

na 

reject  the  extreme  value  distribution  (largest  extreme 
value)  as  the  underlying  distribution.  If  the  sample 
DMAX  value  is  less  than  or  equal  to  DMAX  ,  do  not  reject. 

n  9  ot 


I 


Table  II 

Use  of  Critical  Value  Table 


i 

xi 

F(xi) 

Sn(xi) 

)F(xi)-Sn(xi)  j 

1 

-2.36 

.029 

.067 

.038 

.029 

2 

-1.71 

.074 

.133 

.059 

.007 

3 

-1.23 

.125 

.200 

.075 

.008 

4 

-0.56 

.218 

.267 

.049 

.018 

5 

0.17 

.33  8 

.333 

.00  5 

.071 

6 

1.16 

.506 

.400 

.106 

*  .173  * 

? 

1.54 

.565 

.467 

.098 

.1 65 

8 

1.91 

.619 

.533 

.086 

.152 

9 

2.11 

.  646 

.600 

.046 

.113 

10 

2.63 

.710 

.667 

.043 

.110 

11 

3.37 

•  785 

•  733 

.052 

.118 

12 

3.42 

.789 

.800 

.011 

.056 

13 

3.72 

.814 

.867 

•  053 

.014 

14 

4.23 

.850 

.933 

.083 

.017 

15 

4.92 

.889 

1.00 

.111 

.044 
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VI.  Discussion  of  Results 


The  discussion  of  results  which  follows  is  begun  by- 
covering  the  topic  of  computer  program  validity.  Once 
computer  program  validity  is  established,  the  results  are 
then  presented.  Many  of  the  results  can  be  found  in 
tables  and  figures  in  the  appendices. 

Computer  Program  Validity 

Before  a  major  discussion  of  the  results  begins, 
the  question  of  program  validity  must  be  answered. 

In  order  to  measure  the  quality  of  the  linear  estima 
tion  technique  as  an  estimator  of  the  true  critical  value 
it  is  certainly  necessary  that  the  set  of  critical  values 
used  as  the  standard  set  be  correct.  In  other  words,  the 
validity  of  the  computer  program  generating  the  standard 
technique  critical  values  must  be  verified. 

This  verification  was  initiated  by  generating  in 

f 

the  program  extreme  value  distribution  (largest  extreme 
value)  random  deviates.  The  random  deviates  and  the 
location  and  scale  parameters  were  transformed  to  random 
deviates  and  parameters  appropriate  to  the  two-parameter 
Weibull  distribution  (see  page  25).  Within  the  program, 
the  extreme  value  distribution  (largest  extreme  value) 
cumulative  distribution  function  was  replaced  by  the 
cumulative  distribution  function  for  the  two-parameter 
Weibull,  and  critical  values  were  found  at  n  =  30  for 


significance  levels  .20,  .15,  .10,  .05,  and  .01. 
These  critical  values  should  have  been  two- 


parameter  Weibull  critical  values.  In  order  to  verify 
this  and  thus  the  program  validity,  these  critical 
values  were  compared  to  two-parameter  critical  values 
obtained  at  the  same  sample  size  and  significance  levels 
by  Littell,  McClave,  and  Of fen  (Ref  19:262).  Table  III 
presents  the  two  sets  of  critical  values.  As  Table  III 
indicates,  the  two  sets  of  critical  values  are  almost 
identical,  thus  establishing  the  validity  of  the  compute 
program. 

Selection  of  Plotting  Positions 

In  Appendix  A  is  the  table  of  critical  values  gener 
ated  by  the  standard  technique.  It  was  the  critical 
values  at  n  =  10,  25,  and  40  which  formed  the  standard 
set  of  critical  values,  the  set  of  critical  values  taken 
as  truth. 


Table  III 
Program  Validity 


.20 

.15 

.10 

.05 

.01 

Littell' s 
Results 

.130 

.136 

.  144 

.156 

. 

.179 

Thesis 

Results 

_  .J 

.130 

,137 

.145 

.15? 

.131 
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In  Appendix  B  are  tabled  the  critical  values  generated 
by  the  linear  estimation  technique.  The  critical  values 
in  the  tables  in  Appendix  3  were  formed  by  calculating 
the  DMAX  values  with  the  plotting  position  S  (x)  but 
plotting  the  DMAX  values  with  each  of  the  six  different 
plotting  positions  discussed  in  this  study.  Each  different 
plotting  position  used  to  plot  the  DMAX  values  forms  a 
different  segment  in  the  tables.  Also,  under  each  segment 
is  given  the  mean  absolute  difference  and  the  largest  ab¬ 
solute  difference  related  to  that  set  of  critical  values. 

By  examining  the  various  mean  and  largest  absolute  value 
differences,  it  is  seen  that  the  plotting  positions  defined 
by  Eqs  (15)  and  (17)  at  N  =  750  do  the  best  job  of  reaching 
the  criterion  values  deemed  acceptable. 

It  should  be  noted  that  at  N  =  750  all  the  plotting 
positions  used  to  plot  the  DMAX  values  offer  acceptable 
estimates.  However,  Eqs  (15)  and  (17)  are  equally  better 
than  the  ethers  by  providing  a  mean  absolute  difference  of 
.0012  and  a  largest  absolute  difference  of  .0025.  Since 
Eqs  (15)  and  (17)  provide  equivalent  estimates,  as  measured 
by  mean  absolute  difference  and  largest  absolute  difference, 
the  less  complicated  of  the  two  equations,  Eq  (15),  was 
chosen  as  the  best  plotting  position  with  which  to  plot  the 
DMAX  values. 

At  this  point,  it  nas  been  determined  that  the  linear 
estimation  technique  is  a  "good"  estimation  technique;  how- 
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ever,  substituting  Sn(x)  in  calculating  the  DMAX  values 
with  some  other  plotting  positions  may  improve  the 
estimates . 

In  determining  if  it  would  be  wise  to  replace  Sn(x) 
with  some  other  plotting  position  in  calculating  the 
DMAX  values,  the  plotting  positions  defined  by  Eqs  (15)  > 
( 16) ,  and  (17)  were  selected  as  alternatives.  Thus,  at 
N  =  750,  three  sets  of  critical  values  were  formed  by 
calculating  DMAX  values  with  the  three  alternative 
plotting  positions  and  plotting  all  DMAX  values  with 
Eq  (15).  The  set  formed  by  using  Eq  (15)  in  the  calcu¬ 
lation  of  DMAX  values  produced  a  mean  absolute  differenc 
of  .0114  and  a  largest  absolute  difference  of  .0284. 

The  set  created  by  using  Eq  (16)  had  a  mean  absolute 
difference  of  .0102  and  a  largest  absolute  difference  of 
.0256.  Eq  (17)  formed  a  set  of  critical  values  with  a 
mean  absolute  difference  of  .0032  and  a  largest  absolute 
difference  of  .0080.  These  values  are  definitely  not 
better  or  as  good  as  the  absolute  differences  obtained 
when  Sn(x)  is  used  to  calculate  the  DMAX  values. 

These  results  indicate  that  Eq  (15)  used  to  plot 

the  DMAX  values  and  S  (x)  used  to  calculate  the  DMAX 

n 

values  provide  the  most  effective  linear  estimation 
technique  possible  from  the  plotting  positions  investi¬ 
gated. 


Using  Eq  (15)  and  S  (x)  as  indicated  in  the  previous 
paragraph,  a  table  of  estimated  critical  values  was 
produced  at  N  =  750-  This  table  can  be  examined  in 
Appendix  A. 

The  Power  Study  Results 

In  Appendix  E  can  be  found  tables  indicating  the 
power  of  the  most  effective  linear  estimation  technique 
at  N  =  100,  250,  500,  and  750.  Also,  in  the  same  tables 
are  the  power  values  of  the  standard  technique.  This 
allows  a  direct  comparison  between  the  powers  of  the  stan¬ 
dard  and  linear  estimation  techniques. 

The  power  values  for  the  standard  technique  indicate 
that  the  K-S  test  offers,  among  the  distributions  examined, 
the  best  power  against  the  chi-square  distribution  (1 
degree  of  freedom)  and  the  worst  power  against  the  chi- 
square  distribution  (4  degrees  of  freedom).  As  might  be 
expected,  the  power  against  all  the  distributions  is 
poor  at  n  =  10. 

The  power  values  given  by  the  linear  estimation  tech¬ 
nique  approach  those  of  the  standard  technique  as  N  gets 
larger.  In  order  to  view  the  effect  on  the  power  of  the 
linear  estimation  technique  as  N  increases,  plots  of 
power  versus  N  are  presented  in  Appendix  F.  In  these 
plots,  the  dashed  horizontal  lines  and  the  full  horizontal 
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lines  depict  the  power  associated  with  the  significance 
levels  .01  and  .05,  respectively. 

Examining  the  plots  in  Appendix  F,  it  can  be  seen 
that  each  set  of  plots  related  to  sample  size  n  has 
basically  the  same  form.  This  indicates  that  the  manner 
in  which  the  power  varies  as  N  increases  does  not  depend 
on  the  true  underlying  distribution  of  the  sample  but  on 
the  sample  size  n.  It  is  n  which  determines  if  the 
critical  value  estimated  by  the  linear  estimation  tech¬ 
nique  at  some  N  will  be  smaller  or  larger  than  the  stan¬ 
dard  critical  value.  If  the  critical  value  is  estimated 
to  be  smaller  (larger)  than  the  standard  critical  value, 
the  significance  level  will  be  higher  (lower)  than  the 
standard  and  the  power  will  also  be  higher  (lower)  than 
the  standard. 

"> 

Significance  Levels  * 

Appendix  C  presents  the  proportion  of  rejections 
when  the  hypothesized  distribution,  the  extreme  value 
distribution  (largest  extreme  value),  is  the  actual 
underlying  distribution.  In  other  words,  Monte  Carlo 
generated  significance  levels  are  tabled.  It  can  be 
seen  that  the  standard  technique  does  well  in  meeting  the 
claimed  significance  levels.  The  linear  estimation  tech' 
nique  at  N  =  100  does  poorly;  however,  at  N  =  7p0,  it 


does  well. 


In  Appendix  D  are  plots  of  the  significance  level 
values  versus  N.  The  dashed  and  full  lines  indicate 
the  significance  level  values  for  .01  and  .05,  respec¬ 
tively.  These  plots  present  the  effect  increases  in  N 
have  on  the  sign  'cance  level  values.  Generally,  these 
plots  indicat-3  that  as  N  increases  the  significance  level 
values  approacn  the  standard  significance  levels. 

Chapter  VI  Summary 

This  chapter  began  by  showing  how  the  major  computer 
program  used  in  this  work  was  validated  as  to  giving 
correct  critical  values.  Then,  the  results  of  the  study 
were  presented.  These  results  and  the  conclusions 
stemming  from  them  will  be  presented  in  the  next  and 
final  chapter. 


\ 
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VII.  Conclusions  and  Recommendations 


The  results  introduced  in  the  preceeding  chapter 
leads  to  the  conclusions  offered  here.  In  addition  to 
the  conclusions,  recommendations  for  extensions  of  this 
thesis  are  given. 

Conclusions 

The  following  conclusions  were  drawn  from  the 
results : 

1.  The  standard  set  of  critical  values  that  were 
generated  are  valid. 

2.  The  most  effective  linear  estimation  technique 
is  the  one  which  Eq  (15)  is  used  to  plot  the  DMAX  values, 
Sn(x)  is  used  in  calculating  DMAX,  and  N  =  750. 

3.  The  power  or  significance  levels  are  not  altered 
extensively  by  the  most  effective  linear  estimation  tech¬ 
nique  . 

These  conclusions  can  be  written  in  terms  of  the 
first  two  objectives  set  out  in  Chapter  1. 

1.  The  linear  estimation  technique  does  reduce  the 
DMAX  sample  size  required  to  find  suitable  critical  values. 

2.  Eq  (15)  is  better  than  3n(x)  in  plotting,  the 
DMAX  values;  however,  no  other  plotting  position,  of 
those  studied,  were  better  than  Sn(x)  in  calculating  DMAX. 


45 


And,  of  course,  the  third  and  final  objective  in 
Chapter  1  was  met  by  producing  a  table  of  estimated 
critical  values  using  the  most  effective  linear  estima¬ 
tion  technique. 

Recommendations 

Several  recommendations  for  expanding  this  study 
are  proposed. 

1.  An  investigation  could  be  made  of  the  value  of 
different  interpolation  techniques  to  form  critical 
value  estimates. 

2.  Plotting  positions  not  presented  in  this  thesis 
could  be  examined  for  their  usefulness  in  an  estimating 
technique . 

3.  A  search  could  be  made  to  find  a  functional 
relationship  between  N  and  the  change  in  power  or  signi¬ 
ficance  level. 

4.  An  examination  of  the  linear  estimating  tech¬ 
nique  could  be  made  when  all  the  parameters  of  the  hypo¬ 
thesized  distribution  are  known. 
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Appendix  A 

Tables  for  the  Standard 
Technique  and  the  Most  Effective 
Linear  Estimation  Technique 


Table  V 

Most  Effective  Linear  Estimation  Technique  Critical  Values 


Sample 

Size 

ALPHA  LEVEL 

n 

.20 

•  15 

.10 

.05 

.01 

4 

.  3166 

.3297 

.3419 

.3740 

.4655 

6 

.2699 

.2835 

.2986 

.3234 

•  3714 

8 

.2360 

.2492 

.2685 

.2841 

.3308 

10 

.2205 

.2316 

.2439 

.2607 

.3051 

12 

.1998 

.2090 

.2257 

.2402 

.2747 

14 

.1846 

.1921 

.2037 

.2202 

.2508 

16 

.1736 

.1809 

.1942 

.2101 

.241? 

18 

.1641 

.1713 

.1814 

.1972 

.2398 

20 

.1592 

.  1665 

.17  66 

.1908 

.2228 

25 

.1412 

.1497 

.1584 

.1704 

.1990 

30 

.1322 

.1390 

.1459 

.1570 

.1824 

35 

.1200 

.1251 

.1354 

.1465 

.1733 

40 

.1135 

.1185 

.1275 

.1387 

.1609 

45 

.1065 

.1109 

.1177 

.1274 

.  1447 

50 

.1018 

.1072 

.1135 

.1270 

.140? 

55 

.0963 

.1027 

.1088 

.1172 

.134  2 

60 

.0931 

.0993 

.1048 

.1146 

.1285 
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Appendix  B 


Tables  for  the  Linear  Estimation 
Technique  where  DMAX  Values  are 
Calculated  with  Sn(x) 


Table  VI 


Linear  Estimation  Technique  N  -  100  (a) 


Sample 

Size 

SIGNIFICANCE  LEVEL 

n 

.20 

.15 

.10 

•  05 

.01 

DMAX  plotted  with  S 

(x) 

DMAX  calculated  with 

SnW 

10 

.2220 

.2354 

.248? 

.2576 

.2988 

25 

.142? 

.1521 

.1622 

.1694 

.1907 

40 

.1103 

.1144 

•  12?5 

.1321 

.1445 

LARGEST  DIFFERENCE = 

.0165 

MEAN  DIFFERENCE=  .0049 

DMAX  plotted  with  Eq 

(13) 

DMAX  calculated  with 

Sn(x) 

10 

.2242 

.2357 

.2492 

.2579 

.3072 

25 

.1431 

.1553 

.1639 

.1704 

.1942 

40 

.1118 

.1146 

.127? 

.1326 

.1498 

LARGEST  DIFFERENCE = 

.0112 

MEAN  DIFFERENCE =  .0043 

l  ~  ~  ~  _ 

|  DMAX  plotted  with  Eq 

(14) 

DMAX  calculated  with 

Sn(x) 

10 

.223? 

.2356 

.2491 

.2 578 

.304? 

25 

.1430 

.1545 

.1634 

.1701 

.1932 

40 

.1115 

.1146 

.1276 

.1325 

.1482 

LARGEST  DIFFERENCE^  .0128 


MEAN  DIFFERENCE=  .004 3 


Table  VII 


Linear  Estimation  Technique  N  =  100  (b) 


Sample 

Size 

SIGNIFICANCE  LEVEL 

n 

.20 

.15 

,10 

.05 

.01 

DMAX  plotted  with  Eq 

(15) 

DMAX 

calculated 

with 

sn!x) 

10 

.2234 

•  2356 

.2490 

.2578 

.3030 

25 

.1429 

.15^0 

.1631 

.1699 

.1925 

40 

.1113 

.1146 

.1276 

.1324 

.1472 

LARGEST  DIFFERENCE* 

.0138 

MEAN 

DIFFERENCE* 

.0042 

i  1  "  i 

DMAX  plotted  with  Eq 

(16) 

DMAX 

calculated 

with 

Sn(x> 

10 

.2236 

.2356 

.2491 

.2578 

.3041 

25 

.1430 

.1543 

.1 633 

.1700 

.1929 

40 

.1114 

.1146 

.1276 

.1324 

.1478 

LARGEST  DIFFERENCE- 

.0132 

MEAN 

DIFFERENCE* 

.0043 

l  _ 

DMAX  plotted  with  Eq 

(17) 

DMAX 

calculated  with 

S„(X) 

10 

.2234 

.2356 

.2490 

.2577 

.3030 

25 

.1429 

.1540 

.1631 

.1699 

.1925 

40 

.1112 

.1146 

.1276 

.1323 

.1471 

LARGEST  DIFFERENCE* 

.0139 

MEAN 

DIFFERENCE* 

.0043 

Table  VIII 


Linear  Estimation  Technique  N  =  250  (a) 


Sample 

Size 

SIGNIFICANCE  LEVEL 

n 

..20 

..15  .10 

•.0  5 

-.01 

DMAX  plotted  with  Sn(x)  DMAX  calculated  with  S  (x) 


10 

.2156 

.2267 

.242  5 

.2553 

.2904 

25 

.1340 

. .1460 

.1548 

.1671 

.2015 

40 

.1110 

.1154 

.124? 

.1342 

.1472 

LARGEST 

DIFFERENCE2  .0138 

MEAN 

DIFFERENCE®  .0052 

'  ... 

DMAX  plotted  with  Eq  (13) 

DMAX 

calculated  with 

Sn(x) 

10 

.2159  .2275 

.2427 

.2563 

.3030 

25 

.1403  . 1464 

.1554 

.1678 

.2115 

40 

.1113  .1160 

.1250 

.1353 

.1501 

LARGEST 

DIFFERENCE2  .0112 

MEAN 

DIFFERENCE®  .0040 

1  - — — 

DMAX  plotted  with  Eq  (14) 

DMAX 

calculated  with 

Sn(x) 

10 

.2158  .2271 

.242  6 

.2560 

.3005 

25 

.1402  .1462 

.1553 

•  1675 

.2098 

40 

.1113  .1157 

.1250 

.1347 

.1500 

LARGEST 

DIFFERENCE®  .0110 

MEAN 

DIFFERENCE®  .0042 
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Table  IX 

Linear  Estimation  Technique  N  -  250  (b) 


Sample 

Si  7.  ft 

SIGNIFICANCE 

LEVEL 

n 

.20 

.15 

.10 

.05 

.01 

DMAX  plotted  with  Eq 

(15) 

DMAX 

calculated  with 

Vx> 

10 

.2157 

.2268 

.2426 

.2 557 

.2988 

25 

.1402 

.1461 

.1551 

.1674 

.2087 

40 

.1112 

.1154 

.1249 

.1344 

.1499 

LARGES 1 

?  DIFFERENCE- 

0111 

MEAN 

DIFFERENCE-  .0043 

DMAX  plotted  with  Eq 

(16) 

DMAX 

calculated  with 

sn(x) 

10 

.2158 

.2270 

.2426 

.2559 

.2999 

25 

.1402 

.1462 

.1552 

.1675 

.209^ 

40 

.1113 

.1156 

.1249 

.1346 

.1499 

LARGES 1 

DIFFERENCE-  . 

0111 

MEAN  DIFFERENCE-  .00 4 

Z  . 

DMAX  plotted  with  Eq 

(17) 

DMAX 

calculated  with 

S  n(x) 

10 

.2157 

.2268 

.2426 

.2 557 

.2988 

25 

.1402 

.1461 

.1551 

.1674 

.2087 

40 

.1112 

.1154 

.1249 

.1344 

.1498 

LARGEST  DIFFERENCE-  .0112 


MEAN  DIFFERENCE-  .0043 


Table  X 

Linear  Estimation  Technique  N  -  50°  (a) 


Sample 

Size 

n 


SIGNIFICANCE  LEVEL 


>20 

.15 

.10 

.05 

.01 

|  DMAX  plotted  with 

Sn(x) 

DMAX  calculated  with 

Sn(x)  ! 

10 

.2191 

.2290 

.2425 

.2576 

.3012 

25 

.1413 

.1483 

.1567 

.1674 

.2002  | 

40 

.1119 

.1170 

.1247 

.1368 

j 

1.1546 _ 

LARGEST  DIFFERENCE =  .0064  MEAN  DIFFERENCE^  .0024 


DMAX  plotted  with  Eq  (13)  DMAX  calculated  with  Sn(x) 


10  i 

.2193 

.2291 

.2427 

.2579 

.3072 

25 

.1415 

.1487 

.1568 

.1683 

1 

.2064 

40 

.1121 

.1176 

.1250 

.1386 

.1552 

LARGEST 

DIFFERENCE- 

.0061 

MEAN  DIFFERENCE  =•  .0027 

DMAX  plotted  with  Eq 

(14) 

DMAX  calculated  with  S 

n(x) 

— 

10 

.2192 

.2291 

.24 26 

. 2578 

.3054 

25 

.1415 

.1486 

.1567 

.1680 

.2046 

j  40 

.1121 

.1174 

.1250 

.1381 

.1550 

1 

LARGEST 

1 

DIFFERENCE* 

.0060 

MEAN  DIFFERENCE*  .0025 
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Table  XI 

Linear  Estimation  Technique  N  =  500  (b) 

Sample  i 

Size  SIGNIFICANCE  LEVEL 

i - - 

l 

n 

.20  .15  .10  .05 


DMAX  plotted  with  Eq  (15) 
10  .2192  .2291 

.1414  .1485 

.1120  .1173 


.01 

DMAX  calculated  with  Sn(x) 


.242  6 

.2578 

.3042 

.156? 

.1679 

.2033 

.1249 

.1377 

.1549 

25 

40 


LARGEST  DIFFERENCE-  .0061  MEAN  DIFFERENCE*  .0024 


DMAX 

plotted  with  Eq 

(16) 

DMAX  calculated  with  Sn(x)  | 

10 

.2192 

.2291 

.2426 

.257 8  .3050 

25 

.  1414 

.1486 

.1567 

.1680  .2041 

4, 

t 

t 

.1121 

.1174 

.1249 

i 

.1379  .1550  1 

LARGEST  DIFFERENCE*  . 

0060 

MEAN  DIFFERENCE*  .0025 

DMAX 

plotted  with  Eq 

(17) 

DMAX  calculated  with  Sn(x) 

1 

10  s  .2192 

1 

.2291 

.2426 

.2578 

.3042 

25  .1414 

.1485 

.1567 

.1679 

.2033 

40  .1120 

.1173 

.1249 

.1377 

.1549 

LARGEST  DIFFERENCE-  .0061  MEAN  DIFFERENCE*  .0024 
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Table  XII 

Linear  Estimation  Technique 

N  =  750  (a) 

Sample 

Si  7, e 

SIGNIFICANCE  LEVEL 

n 

.20  .15 

.10 

•  05 

.01 

DMAX  plotted  with  Sn(x) 

DMAX 

calculated 

with  Sn(x) 

10 

.  .  i 

.2202  .2315 

.243? 

.2605 

.3031 

25 

.1412  .1496 

.1583 

.1700 

.1979 

40 

.1134  .1184 

.1275 

.1382 

.1590 

-•-‘'LARGEST  DIFFERENCE*  .002? 

MEAN 

DIFFERENCE; 

-  .0013 

DMAX  plotted  with  Eq  (13) 

DMAX 

calculated 

with  Sn(x) 

10 

.2203  .231? 

.2440 

.2608 

.3061 

25 

.1413  .1498 

.1584 

.1706 

.1996 

40 

.1135  .1185 

.1276 

.1389 

.1626  ; 

LARGEST  DIFFERENCE*  .0035  MEAN  DIFFERENCE*  .0013 


IMAX  plotted  with  Eq  (14)  DMAX  calculated  with  Sn(x) 


10 

.2203 

.2316 

.2439 

.2608 

.3055 

25 

.1412 

.1498 

.1584 

.1705 

.1992 

4o 

.1135 

.1185 

.1275 

.1388 

.1616 

LARGEST  DIFFERENCE*  .0029  MEAN  DIFFERENCE*  .0013 
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Table  XIII 


Linear  Estimation  Technique  N  «  750  (b) 


SIGNIFICANCE  LEVEL 

■ 

.20  -.15 

-.10  -.05 

-.01 

DMAX  plotted  with  Eq  (15) 

DMAX  calculated  with 

Sn(x) 

10 

.2203  .2316 

.2439  .2607 

.3051 

25 

.1412  .1497 

.1584  .1704 

.1990 

40 

.1135  .1185 

.1275  .1387 

.1609 

LARGEST  DIFFERENCE*  .0025 

MEAN  DIFFERENCE*  .0012 

DMAX  plotted  with  Eq  (16) 

DMAX  calculated  with 

Vx> 

10 

.2203  .2316 

.2439  .2607 

.3053 

25 

. 1412  . 1497 

.1584  .1705 

.1991 

40 

.1135  .1185 

.1275  .1387 

.1613 

LARGEST  DIFFERENCE-  .0027 

MEAN  DIFFERENCE*  .0012 

DMAX  plotted  with  Eq  (17) 

DMAX  calculated  with 

s„(x) 

10 

.2203  .2316 

.2439  .2607 

.3050 

25 

.1412  .1497 

.1584  .1704 

.1990 

40 

.1135  .1185 

.1275  .1387 

.1609 

LARGEST  DIFFERENCE*  .0025 

MEAN  DIFFERENCE =  .0012 

Appendix  C 

Table  of  Significance  Levels 
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Table  XIV 


Extreme  Value  (Largest  Extreme  Value) 


Fig  6. 

Extreme  Value  (Largest  Extreme  Value) 


me  Value  (Largest  Extreme 


Appendix  E 

Tables  of  Power  Values 


6? 


Table  XV 
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Table  XVI 


Table  XVII 


XVIII 


Table  XIX 


Appendix  F 

Plots  of  Power  Values 


73 


ma 
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Appendix  G 


Computer  Programs 


-  ♦  „.****-:  • 


DRR,T50Q,CM65000.  TB00627, ROGERS, 44S8 
ATTACH, IHSL, ID-LIBRARY, SN-ASD. 

LIBRARY, INSL. 

FTN5, ANSI-O. 

LSO. 

IEOR 

PROGRAM  THEES2 

C  Hit  NUMBERING  OF  PLOTTING  POSIT  IONS  (PLP)  tUUU 

C  I  I«ITH  ORDERO  RANCOR  DEVIATE  I 

C  I  N-NUHBER  OF  RANDOM  DEVIATES tOR  DNA  STATISTICS)* 

C  I  ilid/N)  t 

C  t  #2i(I-.5)/(N)  » 

C  I  t3i(I-.J)/(N+.4)  I 

C  $  t4»((I/(N+l))+((I-l)/(N-l)))/2  » 

C  »  «5:<I-.375)/(N+.25)  » 

C  I  I6i (I/(N+1))  I 

c  mtmtmmmmmmmmmmmmtmt 

c 

c  tm  condition  codes  Htmmumtmmumtmtum 


C  t  C0ND2: PLOT  D  STATISTICS  USING  PLPI2  t 
C  I  CONDJiPLOT  0  STATISTICS  USING  PLP43  t 
C  t  C0ND4:  PLOT  D  STATISTICS  USING  PLPI4  t 
C  t  C0ND5i  PLOT  D  STATISTICS  USING  PLPI5  t 
C  t  COND61  PLOT  D  STATISTICS  USING  PLPI6  t 


C  I  IF  A  CONDITON  IS  GIVEN  ANY  INTE6ER  VALUE  OTHER  THAN  1, 
C  I  THE  CORRESPONDING  COMPUTATION  HILL  NOT  OCCUR 

C  I  FURTHER,  SNC  MUST  BE  ASSI6NED  AN  INTEGER  VALUE  1 

C  t  FROM  1  TO  5.  THE  VALUE  CORRESPONDS  TO  A  PLP  USED  TO 

C  I  COMPUTE  SN.  1 

c  mmmmmtmmmmmtmmtmmmmstm 

c 

DIMENSION  R (10) , ABDIF (20) , DNA (5000) , SN ( 1 I ) 

INTEGER  REP, SNC 
DOUBLE  PRECISION  DSEED 
DATA  A,B/2I1.0/ 

C-0.0 

NRD-1C 

NAXLIK-SOO 

REP-5000 

SNC-1 

DSEED- 1 234567. ODO 
DO  501  1-1,5000 

C  tm  GENERATES  EITR.  VAL.  RANDOM  DEVIATES  Ml* 

CALL  S6HIB(DSEED,A,NRD,R) 

DO  1  Il-l , NRD 

R(ID— (L06(BIR(I1)+C) ) 

1  CONTINUE 
LOGB— (LQB(B) ) 

THETA-1. 00/A 

c  tmtmmmmutmmmmtmtmmm 

C  Hit  CALCULATES  MAXIMUM  LIKELIHOOD  ESTIMATES  ttlt 

SUM1-0.0 
SUM2-0.0 
SUN3-0.0 
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SUH4-0. 0 
DO  3  I3«ltNRD 
SUH1«R(I3)+SUH1 

3  CONTINUE 

DO  4  14*1 , HAXLIK 
80  5  15*1 f  NRD 
ZX*-(R( 15) /THETA) 

SUN2*R (IS) tEXP (2X) +SUH2 
SUH3*EXP(ZX)+SUH3 

5  CONTINUE 

THETA1- <  SUH1/NRD) - (SUH2/SUH3 ) 

DIF*ABS(THETA1-THETA) 

THETA-THETA1 

IF (DIF. LE. 0. 0000 1 )  60T0100 

4  CONTINUE 

100  00  6  16*1 1 NRD 

ZX1*-(R( 16) /THETA) 

SUM4*EXP ( Z  X 1 ) +SUH4 

6  CONTINUE 

EPS* (-THETA) l(L06(SUH4/NRD)) 

c  ttttttuttttututtttttttmttmmutmttuutt 

C  II!!  GENERATES  DMA  STATISTICS  III! 

CALL  VSRTA(R,NRD) 

SN<1)«0.0 
DO  7  I7*2,NRD+1 
XI7*REAL(I7) 

XNRO*REAL (NRD) 

IF(SNC.EQ.l)  THEN 
SK(I7)*(XI7-1.0)/ (XRRD) 

ENOIF 

7  CONTINUE 

DO  10  J 1*1, NRD 

CDF-EXP4 -EXP ( (EPS-R ( J1 ) ) /THETA) ) 

ABDIF(Jl) *A8S ( CDF-SN  <  J 1 ) ) 

ABDIF ( J 1+NRD) *ABS (CDF-SNI Jl+1 ) ) 

10  CONTINUE 
NRD2*2tNR0 

CALL  V8ATA(A80IF,NRD2) 

DNA(I)*ABDIF(NRD2) 

301  CONTIMJE 

CALL  V8RTA(DNA,REP) 

c  mitttmttttttmmmttmtttimmmm 

C  III!  PRINTS  CRITICAL  VALUES  FOR  5000  REPS  tttl 

PRiNTt,  ’  mmmmmmmttmttmtmtt’ 
PRINT!, 'THE  F0LL0NIN6  CRITICAL  VALUES' 

PRINT!, 'HERE  CALCULATED  FROM  3004’ 

PRINT!, 'DIFFERENT  SETS  OF  ’ , NRD 
PRINT!, 'RANDON  DEVIATES’ 
PRINT!,’!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!’ 
PRINT!, 'THE  80TH1*  ’,DNA(4000) 

PRINT!, ’THE  83TH1*  ’,DNA(4250> 

PRINT!, ’THE  90THX*  ’,DNA(4300) 

PRINT!, 'THE  93THX*  ’,DNA(4750) 

PRINT!,’ THE  99THX*  ’,DNA(4930) 

STOP 
END 
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DRR,  TWO,  CM05000.  T800027, ROGERS, 4458 

ATTACH, 1HSL, ID-LIBRARY , SM-ASD . 

LIBRARY, IRSL. 

FTM5, ANSI-O. 

ISO. 

IEOR 

PROGRAM  CRIVAL 

C  tUt  NUMBERING  OF  PLOTTING  POSITIONS(PLP)  limit 

C  t  I-ITH  OROERO  RANDOM  DEVIATE  I 

C  I  N-NUHBER  OF  RANDOM  DEVIATES (OR  ONA  STATISTICS)! 

C  I  fit (I/N)  $ 

C  »  «2s(I-.5)/(N)  I 

C  »  #3j (I-.3)/(N+.4)  I 

C  II4j((I/(N*1))+((I-I)/(N-1)))/2  I 

C  I  IS)(I-.375)/(N*.25)  t 

C  I  IOiU/IN+I))  I 

c  mmttumtnnmmiimnutmmtnumt 

c 

c  mi  condition  codes  tmmmmimtmmmmmm 


C  I  CONDI i PLOT  D  STATISTICS  USING  PLPI1  I 
C  I  C0ND2IPL0T  D  STATISTICS  USING  PLP«2  I 
C  I  C0ND3i  PLOT  D  STATISTICS  USING  PLPfS  I 
C  I  C0ND4i PLOT  D  STATISTICS  USING  PLPI4  I 
C  I  COM3!  PLOT  0  STATISTICS  USING  PLPfS  I 
C  I  CONDO) PLOT  D  STATISTICS  USING  PLPIO  I 


C  I  IF  A  CONDITON  IS  GIVEN  ANY  INTE6ER  VALUE  OTHER  THAN  1,1 

C  I  THE  CORRESPONDING  COMPUTATION  MILL  NOT  OCCUR  I 

C  I  FURTHER,  SNC  MUST  BE  ASSIGNED  AN  INTEGER  VALUE  ! 

C  I  FROM  1  TO  0.  THE  VALUE  CORRESPONDS  TO  A  PLP  USED  TO  I 

C  I  COMPUTE  SN.  I 

c  mtmmttttttmttitmmmmimttttmmtmm! 

c 

DIMENSION  R (40) , ABDIF (SO) ,DNA(750) , Y (752) , TDNAI7S2) ,SN(41 ) 
INTEGER  REP, DIM, CONDI , CQN02, C0ND3, COM4 , COM3, CONDA, SNC 
DOUBLE  PRECISION  DSEED 
DATA  A.B/211.0/ 

C-0.0 
NRD-40 
INRD-40 
INRD2-INRDI2 
MAILIK-300 
REP-750 
SNC-1 
CONDl-t 
C0ND2-1 
C0ND3-1 
C0ND4-1 
C0M3-1 
CONDO-1 

DSEED- l 234307 ,000 
IIZ-NRB 

IFINRD.BT.S)  GOT0333 
DO  300  NRB-3,IMD 
333  DO  301  I-l.REP 
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DO  2  12*1, INR02 
AB0IF(12>*-1.0 

2  CONTINUE 

DO  SI  LX*1,INRD 
R (LX) *0.0 
31  CONTINUE 

C  III!  GENERATES  EXTR.  VAL.  RANDOM  DEVIATES  till 
CALL  S8NIB(0SEED,A,NRD,R) 

DO  1  11*1, NRS 

R  ( 1 1 )  *-  (L06  (BIR  ( 1 1 )  +0 ) 

1  CONTINUE 
LOSB— (L06<B)> 

THETA*1.00/A 

c  mtmmimmmmmmmtmmtmm 

C  till  CALCULATES  NAXINUN  LIKELIHOOD  ESTIMATES  till 

SUN1*0 
SUM2*0 
SUN3*0 
SUH4*0 

DO  3  13*1, NRD 
SUM1*R(I3)+SUM1 

3  CONTINUE 

DO  4  I4*1,MAXLIK 
DO  3  15-1, NRD 
ZX*-(R( 15) /THETA) 

SUM2*R I 15) IEXP (ZX ) +SUM2 
SUH3*EXP(ZX)+SUM3 

5  CONTINUE 

THETA1*  ( SUH1 /NRD)  - (SUN2/SUH3 ) 

DIF-AB8 (THETAI-THETA) 

THETA-THETA1 

IF(DIF.LE. 0.00001)  60T0100 

4  CONTINUE 

100  DO  &  16*1, NRD 

ZX1«-(R(IA) /THETA) 

SUM4*EXP(ZX1)+SUM4 

6  CONTINUE 

EPS* (-THETA) I (LOO (SUM4/NRD) ) 

c  immmmmmmmmmi 

C  lilt  ORDERS  FIRST  NRD  ELEMENTS  OF  ARRAY  R  till 
DO  77  I0R1*1,NRD-1 
IN»I0R1*1 
DO  7B  I0R2*IN,NRD 
IF  <R< I0R1 ) .LT.R(I0R2) >  60T078 
SN*R(I0R1) 

R(I0R1)*R(I0R2) 

R(I0R2)*SN 
78  CONTINUE 
77  CONTINUE 

c  tmmmmmmmmmmmmmtmmm 

C  tilt  8ENERATES  DMA  STATISTICS  till 

SN<!)«0.0 
DO  7  I7*2,NND*1 
XI7»REAL(I7) 
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XNRD*REAL(NRD) 

IFiSNC.EB. 1)  THEN 
SN(I7)«(XJ7-l,0)/iXNRD) 

ELSE  IFtSNC.E8.2)  THEN 
SN(I7)*( <XI7-1.0>-.5) / <XNRD> 

ELSE  IF (SNC.Efi.3)  THEN 
SM<I7)*<(XI7-1.0>-.3)/ (XNRD+.4) 

ELSE  IFiSNC.E8.4l  THEN 

SN(I7)*(<<XI7-1.0)/(XNRD+1.Q))+I  (XI7-1.0)  /  (XNRD-1 . 0) ) )  /2.0 
ELSE  IFiSNC.EO.S)  THEN 
SN(l7)«<UI7-1.0)-.373)/(XNRD+.25) 

ELSE  IFiSNC.EB. 6)  THEN 
SN<I7)«<XI7-1.0)/iXNRD+1.0l 
ENDIF 

7  CONTINUE 

00  33  J1«1,NR0 

COF-EXP l-EXP( IEPS-R <J1 ) ) /THETA) ) 

ABDIF(J1)«ABS<CDF-SN(J1I) 

ABDIF(J1+NRD)*ABS(CDF-SN(J1+1) ) 

33  CONTINUE 

CALL  VSRTA(ABDIf,INRD2) 

0NA(I)>AB0IF(INR02) 

SOI  CONTINUE 

CALL  VSRTA(DNA,REP) 

c  mmnmmnuttuuumuumtumtm 

C  UU  ESTIMATES  ENDPOINTS.  INTERPOLATES  BETWEEN  Ult 
C  mtm  POINTS,  AND  FINDS  CRITICAL  VALUES  lllllllll 

DIH*REP+2 
DO  8  I8*REP,1,-1 
TDNA(I8+1)*DNA(I8) 

8  CONTINUE 
Y(1)«0.0 
Y(DIH)«1.0 
IFiCONDi.EB.il  THEN 

CALL  PLOTOiY.DIH) 

CALL  ENDP(TONA,r,DIN,REP) 

PRINTI,’  ’ 

CALL  PRNT(REP,NRD,3NC) 

PRINT!,  ’  PLP  ir 
CALL  CRITVA(Y,T0NA,0!H) 

ENDIF 

IF(C0N02.EB.l)  THEN 
CALL  PL0T1 (Y,OIN> 

CALL  ENDP(TONA,Y,OIN,REP) 

PRINTI, ’  ’ 

CAU  PRNT(REP,NRD,3NC) 

PRINT!,’  PLPI2’ 

CALL  CR1TVA(Y,TDNA,0IH) 

ENDIF 

IF(C0ND3.Efl.l>  THEN 
CALL  PL0T2(Y,01N) 

CAU  ENDP(TDNA,Y,DIH,REP) 

PRINTI, ’  ’ 

CALL  PRNT(REP,NRD,SNC) 
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7 


J 


PRINTS,’  PLP#3’ 

CALL  CRITVA(Y,TDNA,DIH) 

ENOIF 

IF (C0N04.EQ. 1 )  THEN 
CALL  PL0T3<Y,DIN) 

CALL  ENOP(TDNA,Y,DIH,REP) 

PRINTS,’  ’ 

CALL  PRNT(REP,NRD,SNC) 

PRINTS,’  PLPI4’ 

CALL  CRITVA(Y,TONA,DIM) 

ENOIF 

IF(CONOS.EO.l)  THEN 
CALL  PLOTS (Y, DIN) 

CALL  EN9P(TDNA,Y,D!N,REP) 

PRINTS,'  ’ 

CALL  PRNT(REP,NRD,SNC) 

PPlHTt,’  PLPI5' 

CALL  CRITVA(Y,TONA,OIH) 

ENOIF 

IF(CQND6.E8.1)  THEN 
CALL  PL0T5(Y,DIM) 

CALL  ENOPITONA,Y,OIN,R£P) 

PRINTS,’  ’ 

CALL  Pm<KP,m,SHC) 

PRINTS,’  PLPI4’ 

CALL  CRITVAtY, TONA,DIN) 

ENOIF 

IF (IZZ.8T.S)  60TQ301 

300  CONTINUE 

301  STOP 
ENO 

c  SSSSSSSSSSSSSSSSSSSSSSSSSSSSSSS 

C  SSSSSSSSS  SUBROUTINES  SSSSSSSSS 

c  sssssssssssssssssssssssssssssss 

C  SSSS  PLOTS  ALQN6  THE  Y-AXIS  USIN6  PLPtl  SSS 

SUBROUTINE  PLOTO(Y,IOIN) 

DIMENSION  Y(IDIH) 

IA-IOIH-1 
DO  9  H0»2, IA 
N0l»H0-1.0 

YtMO)«<H01)/(IDIM-2.0) 

9  CONTINUE 
RETURN 
ENO 

c  sssssssssssssssssssssssssssssssssssssssssss 

c 

C  SSSS  PLOTS  ALOW  THE  Y-AXIS  USIN6  PLPI2  SSS 

SUBROUTINE  P10T1 IV, JOIN) 

DIMENSION  YUDIN) 

IB-IOIff-1 
DO  10  H1*2,!B 
N11«H1-1.0 

Y(Ml)«(NU-.3)/UDIH-2.0) 

10  CONTINUE 
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I 


RETURN 

END 

tttttttttuittttttttttttmtttttttttn 

rn  PLOTS  AL0N6  Y-AXIS  USIN6  PLPI3  Itl 
SUBROUTINE  PL0T2(Y,I0IN) 

DIMENSION  YIIDIN) 

IOIDIH-1 
DO  12  H2*2,IC 
N21*N2-1.0 

Y(H2)*(M21-.3)/UIDIH-2)+.4) 

12  CONTINUE 
RETURN 
END 

tmmmumtmmtummtmut 

tit  PLOTS  AL0N6  Y-AXIS  USING  PLPI4  Ml 
SUBROUTINE  PL0T3(Y,IDIH) 

DIMENSION  YUDIM) 

ID-IDIH-1 
DO  14  N3*2,ID 
H31*M3-1.0 

Y  (M3)  ■<H31/((IDIH-2.0) +1.0)  ♦(M31-1.0)/UIDIM-2.0)-1.0)  1/2.0 
14  CONTINUE 
RETURN 
END 

tttttttttttttntttttttitttitmttttttt 


ttt  PLOTS  ALONG  Y-AXIS  USING  PLP#5  *11 
SUBROUTINE  PL0T4(Y,IDIH) 

DIMENSION  Y(IDIM) 

IE«IDIN-1 
DO  16  N4«2,IE 
M41«M4-1.0 

Y IH4) » (H41-. 375) / U IDIN-2. 0) ♦. 23) 

16  CONTINUE 
RETURN 
END 

tttttttttttttttttttttttttttttttttttttt 

ttt  PLOTS  ALONG  Y-AXIS  USING  PLPI6  It* 
SUBROUTINE  PLQT3<Y,10IN) 

DIMENSION  YUDIN) 

DO  IS  N5»2,1S 

Y (N3)*(M5-1 , 0) / ( ( IDIH-2. 0) +1 . 0) 

IS  CONTINUE 
RETURN 
END 

mtmmmmmmmmmmmt 


It*  ESTIMATES  ENDPOINTS  ttt 
SUBROUTINE  ENOP!X,Y,IV,IN> 
DIMENSION  XUV),YUV> 


4 


SLOPE*  <Y  (2)-Y  (3)  >  /  (X  (2) -JH3)  > 

(2) -SLOPEIX (2) 

Tl*(-B) /SLOPE 
X(1)*T1 

SLOPE*  <Y  < IM) -Y  < IW+1 ) ) / <X  t IM) -X  < IW+1 ) ) 

B1*Y( IK) -SLOPEIX (IK) 

T2*(1.0-BI) /SLOPE 

X(IV)»T2 

RETURN 

END 

c  mmummmttmmu 

c 

C  m  FINDS  CRITICAL  VALUES  III 

SUBROUTINE  CRITVA(Y,TDNA,IDIM) 

DIMENSION  YUDIN), TDNA(IDIH) 

NZA»IDIN-1 
DO  30  J=80,95,5 
DO  31  K*NZAf 1,-1 
IFlYdO.EQ.J/iOO.O)  THEN 
CVAL-TDNA(K) 

PRINT*, 'THE  ’ , J , ’ THI=  ’ , CVAL 
SOTO  30 
ENDIF 

IFIYOO.LT.J/IOO.O)  THEN 
SLOPE* (Y!K)-Y(K+l))/( TDNA <K) -TDNA (K+l ) ) 
Z»YIK)-SLQPEITDNA(K) 

CVAL* ( < J/ 100. 0) -Z ) /SLOPE 
PRINT*, ’THE  ' ,J,’THI*  ’ ,CVAL 
60T030 

ENDIF 

31  CONTINUE 
30  CONTINUE 

DO  32  L*NZA, 1,-1 
IF<Y(L) .EQ.0.99)  THEN 
CVAL«TDNA(L) 

PRINT*, ’THE  99THX*  \CVAL 
SOTO  42 
ENDIF 

IF«Y(L) .LT.0.99)  THEN 

SLQPE«(Y(L)-Y(L+1))/(TDNA(L)-T0NA<LM)) 

Z1«Y(L)-SL0PEITDNA(L) 

CVAL* (0#  99-Z1 ) /SLOPE 
PRINT*,' THE  99THX*  ’ , CVAL 
SOT042 
ENDIF 

32  CONTINUE 
42  RETURN 

END 

c  mitmmmtmmmm 

c 

C  III  PRINTS  HEADINGS  III 

SUBROUTINE  PRNTIJREP,NRD,I8NC) 

print*,  'mum . . 

PRINT*,’  •  OF  REPS*  ’ , JREP 
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PRINTS, ’  »  OF  RAND.  DEV.*  ' , NRD 
PRINTS, ’  SNi  PIPS', ISNC 
RETURN 
END 
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DRR, TSOO, CH63000.  T800627, R06ERS, 4438 
ATTACH,  IHSL,  ID-LIBRARY,  SN-ASD. 

LIBRARY, INSL. 

FTN5, ANSI-O. 

L60. 

4EQR 

PR06RAH  POWER 

c  mtummtmummmmmuttu 

C  I  DIST-lJ  EXT.  VAL.(TYPE1, LARGE  VAL. )  4 

C  t  DIST-2 J  STANDARD  NORMAL  t 

C  I  DIST-Ji  ME I BULL  t 

C  4  DlST-4*  LOG-NORMAL  4 

C  4  DIST-Si  CHI -SQUARED (1  D.F.) 

C  4  DIST-Ai  CHI -SQUARED (4  O.F.) 

C  444444444444444444444444444444444444444 

DIMENSION  RDEVMO) , RANDM1 ) , ABDIF(BO)  ,SNM1) 
DIMENSION  RDEV1 (40) .R0EV2 (40) ,RDEV3(40> 
INTEGER  DIST, REPS, REJ5, REJ 1 
DOUBLE  PRECISION  DSEED.DSEED1 
DOUBLE  PRECISION  DSEES2, DSEED3 
DATA  A, B,C,D, E/241. 0,240.0, 1.0/ 

DIST-3 

IREPS-230 

PRINT*,'  REPS*’ , IREPS 
PRINT4,'  ♦♦  DISTRIBUTION  *’,DIST,' 

DSEED* 1 234567. ODO 

DSEE01-4567123.0D0 

DSEED2-3451267.0D0 

DSEED3-23476S1 . ODO 

NRD-40 

NRS2-NRD42 

REJ3-0 

REJ 1*0 

DO  98  REPS-1,2500 
C  444444144444444444444 

C  4  GENERATES  DIST. II  4 

C  144444444444444444444 

IF (DIST. EG. 1)  THEN 
CALL  66NIB(DSEED,A,NRD,RDEV) 

D01  Il-l.NRD 

RDEV  ( 1 1 ) — (LOG  (B4RDEV  ( 1 1 )  +0 ) 

1  CONTINUE 
ENDIF 

C  444444444444444444444 

C  4  GENERATES  DIST. 12  4 

C  1444(1444444444444441 

IF (DIST. EB. 2)  THEN 
CALL  S6NML(08EED,NR0,R0EV) 

ENDIF 

C  (44444444444444(44444 

C  4  GENERATES  DIST.I3  4 

C  444444444441114444444 

IF (SIST.Efl.3)  THEN 
CALL  86NII(D8EED,A,NRD,RDEV) 


END  IF 

c  mmmttmttttttt 

C  I  GENERATES  BIST. 14  « 

c  mmtmtmmmt 

IFIDIST.ED.4)  THEN 
CALL  GGNLSIDSEED, NRD,D, £,  RAND) 

DOW  J»1,NRD 
RDEV(J)*RAND(J) 

99  CONTINUE 
ENDIF 

c  mtmmtmtmm 

C  t  GENERATES  DISTtS  4 

c  mmtmmtmm 

IF(DIST.ES.S)  THEN 
CALL  S6NNL ( DSEED, NRD,  ROE V) 

D022  I 22*1, NRD 

RDEV  ( 122) »RDEV  (122)  IRDEV  ( 122) 

22  CONTINUE 
ENDIF 

c  mmtmmtmm 

C  t  GENERATES  DISTI6  I 

c  mttmtttttttttttt 

IF(DIST.Efl.G)  THEN 
CALL  66NHL ( DSEED , NRD , RDEV ) 

CALL  68NNL ( DSEED 1 , NRD , RDEV 1 ) 

CALL  S6NHL (DSEED2,NRD,RDEV2) 

CALL  66NHL ( DSEED3, NRD,  R0EV3 ) 

009  19*1, NRD 

RDEV<I9)*R0EV<I9)tt2+R0EVl<I9)tt2+RDEV2(I9)tt2+RDEV3(I9)tt2 
9  CONTINUE 
ENDIF 

c  mttttttttttttttttttttttttttttttt 

C  I  CALCULATES  PARANETER  ESTIMATES  t 

c  ttmmmmtttmtmmtmm 

EPS*0.0 
THETA-1.0 
SUM 1*0.0 
3UN2-0.0 
SUN3-0.0 
3UH4-0.0 
003  13*1, NRD 
SUN1*RD£V(I3)+8UN1 

3  CONTINUE 
D04  14*1,1000 

DOS  13*1,  NRD 
21*- IRDEV (I 5) /THETA) 

SUM2*RDEV ( IS) tEJ P ( ZZ ) +SUM2 
SUH3«E3PIZI)*SUH3 
S  CONTINUE 

THETA1«(SUH1/NRD)-(SUN2/3UN3> 

DIF-AB8 (THETA1-THETA) 

THETA*THETAl 

IF < DIF .LE. 0.00001 )  GOTO  100 

4  CONTINUE 
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100  000  16*1 , NRO 

m*-<RDEVII6) /THETA) 

SUM4*EXP  <  2X 1 ) +SUM4 

6  CONTINUE 

EPS* ( -THETA) $ ( LOS ( SUH4/NR0 ) ) 

c  mttmumttmmmntmmmmmtm 

C  t  ORDERS  RDEV,  FINDS  ONAI,  COUNTS  REJECTIONS  t 

c  . . . 

CALL  VSRTA(RDEV,NRD> 

SN(1)*0.0 
D07  I7-2,NRD+1 
XI7*R£AL1 17) 

XNRO-REAL(NRD) 

SN(I7)*(XI7-1.0) / (XNRD) 

7  CONTINUE 
0030  10*1, NRO 

CDF*EXP(-EXP( (EPS-RDEV(ID) ) /THETA) ) 

ASDIF ( ID) >ABS (CDF-SN ( ID) ) 

ABDIF ( ID+NRD) >ABS (CDF-SN ( ID* 1 » 

SO  CONTINUE 
CALL  VSRTA(ABDIF,NRD2) 

DHAX*ABDIF(NRD2) 

IF(DNAX.6T. 0.1344)  REJ3-REJ5+1 
JF(DHAI.6T. 0.1499)  REJ1«REJ1+1 
98  CONTINUE 

IREJ3*R£AL(REJ3) 

XREJ1*REAL<REJ1) 

PERC3*XRE J3/2S00. 0 
PERCl*XREJl/2300.0 

print*,  •mtttmumtv 

PRINT*,*  NRD*  *,NRD 
PRINT*, ***************** 

PRINT*,’  REJ3-  *,REJ5,*  REJ1*  *,REJ1 
PRINT*,’ *********** 

PRINT*,’  -SN1-  ’ 

PRINT*,’  ALPHA*. OS’ 

PRINT*, ’I  REJECTED*  ’,PERC5 
PRINT*,’  ' 

PRINT*,’  ALPHA*. 01’ 

PRINT*, ’*  REJECTED  «’,PERCI 
PRINT*,’  ’ 

PRINT*, ’**(***(***’ 

STOP 

END 
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ORR. T400, CH65000.  T800627, ROGERS. 4438 
ATTACH, ML, IO*LIBRARY,SN*ASD. 

LIBRARY, IHSL. 

FTN5,ANSI*0. 

LBO. 

IEOR 

PROGRAM  POWER A 

c  mtmmtuutmtmummutmu 

C  t  DIST-li  EXT.  VAL. (TYPE l, LARGE  VAL. )  I 


C  t  0IST>2>  STANDARD  NORMAL  I 
C  I  DIST*3i  WE I  BULL  4 
C  I  DISTWi  LOG-NORMAL  t 
C  t  DIST*5i  CHI-SOUAREDI1  D.F.)  I 
C  I  D1ST*6»  CHI 'SQUARED (4  D.F.)  I 


c  tmtmtuttmmmmmtmmtm 

DIMENSION  RDEV (40) , RAND ( 41 ) , ABDIF3 (80) , ABDIFS <B0) , ABDIF6 (80) 
DIMENSION  SN3(41),SN5(41),SN6(41) 

DIMENSION  RDEV1 (40),RDEV2(40),RDEV3(40) 

INTEGER  BIST, REPS, REJ33,REJ31,REJ33,REJ51,REJ65,REJ61 
DOUBLE  PRECISION  DSEED,DSEEDl 
DOUBLE  PRECISION  DSEED2, DSEED3 
DATA  A, B,C,D, E/211. 0,210.0, 1.0/ 

0IST«1 

PRINT*, ’  ♦♦  DISTRIBUTION  »’,DIST,’ 

DSEED* t 234567. ODO 

DSEED 1 *45671 23. ODO 

DSEED2*343 1267.  ODO 

DSEED3*2347631.0D0 

NRD*40 

NRD2«NRDI2 

REJ33*0 

REJ33*0 

REJ31-0 

REJ31-0 

REJ65-0 

REJ61*0 

DO  98  REPS* 1,2300 

c  mmtmuttmmt 

C  I  GENERATES  DIST.I1  t 

c  tmmmitmmttt 

IF (DIST.EQ, 1 )  THEN 
CALL  66WIB10SEED,A,NRD,RDEV) 

001  11*1, NRD 

RDEV(I1)*-(L08(BIRDEV(I1)+C)) 

1  CONTINUE 
END  IF 

c  tttmtttttmmtm 

C  I  GENERATES  DIST.I2  I 

c  mitmittmtmm 

IF (BIST, £8.2)  THEN 
CALL  86NML ( DSEED , NRD , RDEV ) 

ENDIF 
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c  mumtmtmmu 

C  t  GENERATES  DIST.I3  I 

c  mtmmttttuutu 

IFlDIST.EB.3)  THEN 
CALL  S6N I 8 ( DSEED , A , NRD , ROE V ) 

ENDIF 

c  mttmunmmm 

C  I  GENERATES  BIST. 14  I 

c  mmtumumtm 

IFIDI8T.ES.4)  THEN 
CALL  GGNLGlOSEED, NRD, D,E, RAND) 

0099  J*1,NR0 
RDEV(J)«RANO(J) 

99  CONTINUE 

ENDIF 

c  tmtttmmtmmm 

C  I  GENERATES  DIST.iS  t 

c  immmmmttm 

IF (DIST.ES.  3)  THEN 
CALL  86MHL <DSEEDyNRD,ROEV) 

D022  122*1, NRD 

RDEV ( 122) »R0EV ( 122) tROEV (122) 

22  CONTINUE 

ENDIF 

c  mumtmuimm 

C  *  GENERATES  DIST.I6  » 

c  mmmmumtm 

IF (DIST.ES. A)  THEN 
CALL  66NHL(DSEED, NRD, RDEV) 

CALL  GGNRL ( DSEED 1 , NRD , R0EV1 ) 

CALL  GGNNL ( DSEED2, NRD, RDEV2) 

CALL  66MH.  ( DSEED3,  NRD,  RDEV3 ) 

D09  19*1, NRD 

RDEV  ( I9)*RDEV(  19)  U2+RDEV1  ( 19)  M2+R0EV2  (19)  U2+RDEV3U9)  U2 
9  CONTINUE 
ENDIF 

c  ttmmttmtmtmtttmtuuu 

C  4  CALCULATES  PARAMETER  ESTIMATES  I 

c  mtmtmtmuummmmtu 

EPS-0.0 
THETA* 1.0 
SUM1*0.0 
SUM2*0.0 
SUM3*0.0 
SUM4-0.0 
003  I 3*1, NRD 
SUM1*RDEV(I3)+SUM1 
3  CONTINUE 
004  14*1,1000 
DOS  13*1, NRD 
ZX*-(RD£V(I5) /THETA) 

SUM2-RDEV (IS) tEXP(ZX) *SUM2 
SUM3*EXP(IX)+SUM3 

3  CONTINUE  1Uji 


« 


THETA1* (SUN1/NR0) - (SUH2/SUH3) 

DIF«ABS ( THETA 1 -THETA ) 

THETA*THETA1 

IF (DIF. LE. 0. OOOOl )  SOTO  100 
4  CONTINUE 
100  006  16*1,HRD 

m»-(ROEVUA)  /THETA) 

SUN4*£XP(ZX1)+SUH4 
A  CONTINUE 

EPS-  (-THETA)  t  (L06  (SUR4/NR0)  > 

c  tmmumtmmutmttmmmtmtmtt 

C  t  ORDERS  RDEV,  FINDS  DHAX,  COUNTS  REJECTIONS  I 

c  mttmmmmmmmmmmmttmm 

CALL  VSRTA(RD£V,NRO) 

SN3ID-0.0 
SN5(11*0.0 
SN6(1)*0.0 
D07  I7*2,NRD+1 
XI7-REAL ( 17) 

XNRD-REAL(NRD) 

SN3 (I7)*( 1117*1.0) -.3)/ (INRD+, 4) 
SN3(I7)*((XI7-1.0)‘.375)/UNRD+.23> 

SNA( 17)- (XI7-1 . 0) / ( XNRD+ 1.0) 

7  CONTINUE 
D030  ID«1,NRD 

CDF-EXP (-EXP( (EPS-RDEVI ID) ) /THETA) ) 

AB0IF3 ( ID) *ABS (CDF-SN3 ( ID) ) 

ABDIF3 ( ID+NRD) -ABS (CDF-SN3 ( ID+1 ) ) 

A8DIF3 (ID) -ABS (CDF-SN3 1 ID) ) 

AB0IF5< ID+NRD) *ABS(CDF-SN3 ( ID+1 ) ) 

ABDIFA ( ID) -ABS (CDF-SN6 ( ID) ) 
ABDIF6(ID)*ABS(CDF-SN6(ID+1)) 

SO  CONTINUE 

CALL  VSRT A  <  ABD I F3, NRD2 ) 

CALL  VSRTA(ABDIF3|NRD2) 

CALL  VSRTA<A80IF6,NRD2) 

DHAX3-ABDIF3 (NRB2) 

DHAX  3- ABD I F3 ( NRD2 ) 

DHAX6»ABDIF6(NRD2> 

IF  (DIWI3.  ST.  0. 1 407)  REJ33-REJ33+1 
IF  (DIWI3.  ST.  0.1659)  REJ31-REJ3U1 
IF (DHAX5.8T.0. 1403)  REJ53«REJ33+1 
IF (DNAX5. ST. 0.1637)  REJ31*REJ31H 
IFI0NAX6. ST. 0.1404)  REJ63-REJ65+1 
IF (DHAX6.8T. 0.1637)  REJ61*REJ61M 
98  CONTINUE 

XREJ33-REALIREJ35) 

XREJ31-REAL (REJ31 ) 

IR£J55»R£AL(R£J53) 

IREJSl-REALtREJSl) 

XREJ6S-REAL (REJ65) 

XREJ6 1 -REAL ( RE J6 1 ) 

print*/ mttmtmtm' 

P33-XREJ3S/2300. 0 
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P31*XREJ31/2500.0 
P35MREJ55/2500.0 
P51*XR£J51/2500,0 
P65*XREJA3/2300.0 
PAWREJ61/2500.0 
PRINT!,’  NRD*  ’,NRD 

print*,  'tmmtmuM' 

PRINT!,’  ’ 
PRINT!,’!!!!!!!!!!’ 
PRINT!,’  -SN3-  ’ 

PRINT!,’  ALPHA-.05* 
PRINT!,’  Z  REJECTED*  ’,P33 
PRINT!,’  ’ 

PRINT!,’  ALPHA*. 01’ 
PRINT!,’  X  REJECTED*  ’,P31 
PRINT!,’  ’ 

PRINT!, ’I!!!!!!!!!’ 
PRINT!,’  -SN3-’ 

PRINT!,’  ALPHA*. 03’ 
PRINT!,’  I  REJECTED*  ’,P33 
PRINT!,’  ’ 

PRINT!,’  ALPHA*. Of 
PRINT!,’  Z  REJECTED*  \P51 
PRINT!,’  ’ 
PRINT!,’!!!!!!!!!!' 

PRINT!,’  -SNA-’ 

PRINT!,’  ALPHA*. 03’ 

PRINT!,’  Z  REJECTED*  ’,PA3 
PRINT!,’  ’ 

PRINT!,’  ALPHA*. Of 
PRINT!,’  Z  REJECTED*  ’,PA1 
STOP 
END 
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